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Abstract 


Adaptive  processing  of  RADARSAT-2  synthetic  aperture  radar  (SAR)  imagery  of  known  ships  ac¬ 
quired  near  the  Strait  of  Gibraltar-  provides  estimates  of  the  ship  radial  speed  along  the  radar  line  of 
site  and  the  ship  aspect  angle,  which  can  be  combined  to  provide  an  estimate  of  the  ship  velocity  in 
the  absence  of  a  visible  ship  wake.  Using  a  Doppler  offset  algorithm,  it  was  found  that  the  standard 
deviation  of  the  radial  speed  estimates  was  0.85  m/s  for  26  ships  in  two  HH  scenes,  compared  to 
0.3  m/s  using  along-track  interferometry  (ATI)  algorithms,  and  was  1.3  m/s  for  9  ships  in  a  VV 
scene,  compared  to  0.6  m/s  using  ATI  algorithms.  It  was  shown  that  ship  aspect  angle  estimation 
can  be  complicated  by  wave-induced  ship  motion,  which  can  smear  or  change  the  orientation  of  the 
imaged  ship  signature.  Validation  was  carried  out  against  coastal-received  Automatic  Identification 
System  (AIS)  data.  It  is  apparent  that  ship  velocity  estimation  is  feasible  for  single  aperture  SARs 
in  the  absence  of  a  visible  ship  wake  signature,  which  could  be  useful  for  tracking  ships  that  are  not 
transmitting  AIS  data. 


Resume 


Le  traitement  adaptatif  d'images  de  navires  connus  obtenues  par  un  radar  a  synthese  d'ouverture 
(SAR)  de  RADARSAT-2  pres  du  detroit  de  Gibraltar-  donne  des  estimations  de  la  vitesse  radiale 
de  navires  le  long  de  l'axe  de  visibilite  directe  du  radar  et  de  Tangle  d'aspect  des  navires,  lesquelles 
peuvent  etre  combinees  pour  obtenir  une  estimation  de  la  vitesse  des  navires  en  Tabsence  d'un 
sillage  visible  des  navires.  A  l'aide  d'un  algorithme  de  decalage  Doppler,  on  a  trouve  que  l’ecart- 
type  des  estimations  de  la  vitesse  radiale  etait  de  0,85  m/s  dans  le  cas  de  26  navires  dans  deux 
scenes  HH,  par  rapport  a  un  ecart-type  de  0,3  m/s  obtenu  a  l'aide  d'algorithmes  d'interferometrie 
longitudinale,  et  de  1,3  m/s  dans  le  cas  de  9  navires  en  mode  VV,  par  rapport  a  0,6  m/s  obtenu  avec 
l'algorithme  d'interferometrie.  On  a  montre  qu'une  estimation  de  Tangle  d'aspect  d'un  navire  peut 
etre  compliquee  par  le  mouvement  du  navire  induit  par  les  vagues,  qui  peut  brouiller  l'image  ou 
modifier  Torientation  de  la  signature  du  navire  qui  parait  a  l'image.  Une  validation  a  ete  effectuee 
par  rapport  a  des  donnees  du  systeme  d'identification  automatique  (SIA)  rccucs  d'une  station  cohere. 
De  toute  evidence,  il  est  possible  d'estimer  la  vitesse  d'un  navire  pour  les  radars  SAR  a  une  seule 
ouverture  en  Tabsence  de  la  signature  d'un  sillage  visible  du  navire,  qui  pourrait  s'averer  utile  a  la 
poursuite  des  navires  qui  n'emettent  pas  de  donnees  SIA. 
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Executive  summary 


Estimation  of  ship  velocity  by  adaptive  processing  of  single 
aperture  RADARSAT-2  data 

Paris  W.  Vachon,  Marina  V.  Dragosevic;  DRDC  Ottawa  TM  2009-109;  Defence  R&D 
Canada  -  Ottawa;  August  2009. 

Background:  Adaptive  processing  of  synthetic  aperture  radar  (SAR)  imagery  of  known  ships  us¬ 
ing  the  CHip-based  Adaptive  SAR  Processor  (CHASP)  developed  at  DRDC  Ottawa,  is  tested  and 
demonstrated  on  three  RADARSAT-2  data  sets  that  were  acquired  near  the  Strait  of  Gibraltar  in 
Moving  Object  Detection  Experiment  (MODEX)  mode.  The  data  were  acquired  in  this  mode  so 
that  single  beam  range-compressed,  azimuth-uncompressed  Phase  History  Product  (PHP)  with  suf¬ 
ficient  fidelity  could  be  created  by  using  beam-forming,  and  so  that  estimated  single  beam  ship 
velocities  could  be  compared  with  estimates  from  dual-aperture  SAR  data.  The  analysis  focuses 
on  estimation  of  the  ship  radial  speed  along  the  radar  line  of  site  and  the  ship  aspect  angle.  These 
two  values  can  be  combined  to  provide  an  estimate  of  the  ship  velocity  in  the  absence  of  a  visible 
ship  wake.  The  radial  speed  is  estimated  by  a  Doppler  offset  algorithm  that  was  previously  devel¬ 
oped  and  demonstrated  for  RADARS  AT- 1  Fine  mode  data.  A  well-calibrated  physical  model  of  the 
Doppler  centroid  (DC)  serves  to  reduce  the  radial  speed  bias  in  the  presence  of  surface  currents. 
The  ship  aspect  angle  may  be  estimated  from  the  ship  signature  in  the  image  domain  using  sev¬ 
eral  algorithms,  including  region  growing  and  an  Eigen-value  decomposition.  The  ships  considered 
were  validated  using  coastal-received  Automatic  Identification  System  (AIS)  data. 

Results:  The  data  sets  include  two  HH  and  one  VV  acquisitions.  It  was  found  that  the  standard 
deviation  of  the  radial  speed  estimates  was  0.85  m/s  for  26  ships  in  the  HH  scenes,  compared  to 
0.3  m/s  that  was  achieved  using  along  track  interferometry  (ATI)  algorithms  on  the  MODEX  mode 
data.  The  performance  was  worse  for  the  VV  scene,  achieving  1.3  m/s  for  9  ships,  compared  to 
0.6  m/s  achieved  with  ATI  algorithms.  A  calibrated  physical  model  for  the  Doppler  centroid  is 
essential  for  the  single-beam  technique  to  be  successful. 

In  this  work,  a  customized  single  look  complex  (SLC)  product  was  generated  from  the  MODEX 
mode  data.  The  analysis  shows  that  SLC  could  be  a  good  starting  point  for  this  analysis,  but  a  special 
product  with  a  large  processed  bandwidth  and  no  azimuth  look  weighting  would  be  required. 

Ship  aspect  angle  estimation  can  be  complicated  by  wave-induced  ship  motion,  which  will  smear 
or  change  the  orientation  of  the  imaged  ship  signature. 

Significance:  Ship  velocity  estimation  is  feasible  for  single  aperture  SARs  in  the  absence  of  a  visi¬ 
ble  ship  wake  signature.  From  an  operational  perspective,  this  would  require  access  to  a  special  SLC 
product.  Derived  velocities  could  be  used  to  forecast  the  location  of  ships  that  arc  not  transmitting 
AIS  data. 


DRDC  Ottawa  TM  2009-109 


Sommaire 


Estimation  of  ship  velocity  by  adaptive  processing  of  single 
aperture  RADARSAT-2  data 

Paris  W.  Vachon,  Marina  V.  Dragosevic  ;  DRDC  Ottawa  TM  2009-109  ;  R  &  D  pour  la 

defense  Canada  -  Ottawa ;  aout  2009. 

Introduction  :  Le  traitement  adaptatif  d'images  de  navires  connus  par  un  radar  a  synthese  d'ouverture 
(SAR)  a  l'aide  du  processeur  SAR  adaptatif  a  puce  (CHASP),  mis  au  point  par  RDDC  Ottawa, 
a  ete  mis  a  l'essai  et  a  fait  l'objet  d'une  demonstration  faisant  appel  a  trois  jeux  de  donnees  de 
RADARSAT-2  acquises  pres  du  detroit  de  Gibraltar-  en  mode  d'experience  de  detection  d'objets 
mobiles  (MODEX).  Les  donnees  ont  ete  acquises  dans  ce  mode  pour  qu'il  soit  possible  de  creer 
assez  fidelement  un  produit  de  discordance  de  phase  (PHP)  a  compression  de  la  portee  d'un  seul 
faisceau  et  sans  compression  de  l'azimut,  au  moyen  de  la  mise  en  forme  des  faisceaux,  pour  qu'il 
soit  possible  de  comparer  les  vitesses  estimees  de  navires  a  un  seul  faisceau  a  des  estimations 
faites  a  partir  de  donnees  de  radars  SAR  a  deux  ouvertures.  L'analyse  porte  principalement  sur 
l'estimation  de  la  vitesse  radiale  de  navires  le  long  de  l'axe  de  visibility  directe  du  radar  et  de  Tangle 
d'aspect  des  navires.  II  est  possible  de  combiner  ces  deux  valeurs  pour  donner  une  estimation  de 
la  vitesse  des  navires  en  l'absence  d'un  sillage  visible  des  navires.  La  vitesse  radiale  est  estimee  a 
l'aide  d'un  algorithme  du  decalage  Doppler  qui  a  ete  elabore  et  a  fait  l'objet  d'une  demonstration 
pour  les  donnees  de  RADARSAT-1  en  mode  de  resolution  fine.  Un  modele  materiel  bien  etalonne 
du  centroide  Doppler  sert  a  la  reduction  de  Terreur  systematique  de  la  vitesse  radiale  en  presence  de 
courants  de  surface.  L'angle  d'aspect  des  navires  peut  etre  estime  a  partir  de  la  signature  des  navires 
dans  la  representation  en  image  a  l'aide  de  plusieurs  algorithmes,  y  compris  la  fusion  de  regions  et 
une  decomposition  de  valeurs  propres.  Les  navires  en  question  ont  ete  valides  a  l'aide  de  donnees 
du  systeme  d'identification  automatique  (SIA)  rccucs  d'une  station  cotiere. 

Resultats  :  Les  jeux  de  donnees  comprennent  des  donnees  obtenues  en  deux  scenes  d'acquisition 
HH  et  un  scene  d'acquisition  VV.  On  a  trouve  que  Tecart-type  des  estimations  de  la  vitesse  radiale 
etait  de  0,85  m/s  dans  le  cas  de  26  navires  en  mode  HH,  par  rapport  a  un  ecart-type  de  0,3  m/s  obtenu 
a  l'aide  d'algorithmes  d'interferometrie  longitudinale  appliques  aux  donnees  en  mode  MODEX.  Le 
rendement  etait  pire  en  mode  VV,  l'ecart-type  etant  de  1,3  m/s  dans  le  cas  de  9  navires,  par  rapport 
a  0,6  m/s  obtenu  avec  Talgorithme  d'interferometrie.  Un  modele  materiel  etalonne  pour  le  centroide 
Doppler  est  essentiel  au  succes  de  la  technique  faisant  appel  a  un  seul  faisceau. 

Dans  le  present  travail,  un  produit  en  format  SLC  (single  look  complex)  adapte  a  ete  genere  a 
partir  des  donnees  en  mode  MODEX.  L'analyse  montre  que  le  format  SLC  pourrait  constituer  un 
bon  point  de  depart  de  l'analyse,  mais  il  faudrait  un  produit  particulier  ayant  une  grande  largeur  de 
bande  de  traitement  et  sans  ponderation  de  Tobservation  de  l'azimut. 

L'estimation  de  Tangle  d'aspect  de  navires  peut  etre  compliquee  par  le  mouvement  des  navires  induit 
par  les  vagues,  qui  brouille  l'image  ou  modifie  Torientation  de  la  signature  des  navires  qui  parait  a 
l'image. 


IV 
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Portee  :  II  est  possible  d'estimer  la  vitesse  d'un  navire  pour  les  radars  SAR  a  une  seule  ouverture  en 
l'absence  de  la  signature  d'un  sillage  visible  du  navire.  Du  point  de  vue  operationnel,  cela  exigerait 
Faeces  a  un  produit  particulier  en  format  SLC.  Les  vitesses  derivees  pourraient  servir  a  la  prevision 
de  l'emplacement  de  navires  qui  n'emettent  pas  de  donnee  SIA. 
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1  Introduction 


In  earlier  work,  the  CHip-based  Adaptive  SAR  Processor  (CHASP)  [1]  was  used  to  adaptively 
process  RADARSAT-1  data,  acquired  in  the  Fine  Beam  mode,  and  to  estimate  significant  parameters 
of  the  detected  ships  [2],  Some  of  the  estimation  methods  arc  now  being  tested  on  RADARS AT-2 
data  with  similar  properties. 

This  work  is  primarily  focused  on  the  estimation  of  two  parameters,  namely  the  radial  speed  and  the 
ship  aspect  angle.  Radial  speed  is  defined  as  the  projection  of  the  ship’s  velocity  onto  the  radar  line 
of  site  (LOS).  Its  positive  direction  is  away  from  SAR.  The  aspect  angle  is  defined  as  the  relative 
heading  of  the  ship,  measured  in  the  ground  plane  stalling  from  the  ground  projection  of  the  LOS. 
Its  positive  direction  is  the  same  as  for  measuring  squint,  i.e.,  rotation  about  a  downwards  axis  is 
taken  to  be  positive. 

Radial  speed  and  relative  heading  uniquely  define  the  velocity  of  the  ship  since  the  acquisition 
geometry  is  known.  However,  the  estimation  accuracy  depends  on  many  factors.  In  this  work,  radial 
speed  is  estimated  by  the  same  Doppler  offset  algorithm  previously  developed  for  RADARSAT-1. 
It  has  been  shown  that  a  well-calibrated  physical  model  of  the  Doppler  centroid  (DC)  plays  an 
important  role  in  reducing  the  bias  of  the  radial  speed  estimates  in  the  presence  of  surface  currents 
[3],  Therefore,  special  attention  is  paid  to  the  DC  model  for  RADARS  AT-2.  The  ship  aspect  angle 
is  estimated  by  the  region  growing  algorithm  applied  to  the  image  signature,  which  is  a  very  simple 
technique  initially  developed  to  provide  rough  estimates  of  the  vessel’s  size  and  orientation.  A  new 
algorithm,  based  upon  Eigen  value  decomposition  of  ship  signature,  is  developed  and  also  used  in 
this  work.  Estimation  errors  of  these  algorithms  arc  analysed  in  more  detail  than  in  previous  work. 

All  SAR  data  used  in  this  work  have  been  acquired  in  the  special  split-antenna  (dual  channel) 
mode,  but  have  been  converted  to  single  aperture  (single  channel)  SAR  data  by  beamforming.  The 
main  reason  for  this  approach  is  that  only  dual  channel  RADARSAT-2  data  are  available  in  their 
raw  (unfocused)  form.  As  a  result  of  this  approach,  it  is  also  possible  to  compare  the  accuracy 
of  the  applied  single-aperture  adaptive  algorithms  to  the  accuracy  of  the  well-known  dual  aperture 
techniques  [4],  [5], 

While  there  is  a  huge  body  of  open  literature  on  velocity  estimation  with  dual  channel  SAR,  velocity 
estimation  with  single  aperture  SAR  remains  a  challenging  topic  with  very  few  real-world  solutions. 
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2  Properties  of  RADARSAT-2  SAR  data  collection 


The  SAR  data  collection  used  in  this  work  is  similar  to  the  RADARSAT-1  Fine  Beam  Mode  data 
collection  that  was  the  basis  of  the  previous  work.  The  wavelength,  incidence  angle  and  resolution 
arc  all  very  similar.  However,  RADARSAT-2  has  some  different  properties  related  to  its  steered 
phased-array  antenna.  Steering  is  used  to  create  close-to-zero  DC  conditions  for  the  acquired  signal. 
Also,  there  is  the  capability  of  using  the  antenna  in  the  split  mode,  which  is  the  key  feature  exploited 
in  the  Moving  Object  Detection  Experiment  (MODEX)  [6],  Unfortunately,  the  RADARSAT-2  data 
policy  does  not  make  the  raw  signal  data  available  for  any  other  use,  even  though  many  applica¬ 
tions  could  benefit  from  the  Phase  History  Product  (PHP),  that  is,  range  compressed  and  azimuth 
uncompressed.  Therefore,  PHP  used  in  this  work  is  generated  by  beamforming  of  the  two-channel 
PHP  produced  in  MODEX  mode. 

2.1  Doppler  centroid  model  for  RADARSAT-2 

Under  the  assumption  that  attitude  angles  (yaw,  (|>y,  pitch,  typ,  and  roll,  <{>,.)  are  small,  a  DC  expression 
can  be  derived  in  the  closed  form.  Linearization  with  respect  to  yaw  and  pitch  is  justified  because 
these  angles  arc  small.  In  the  first  order  approximation,  roll  can  be  absorbed  into  the  off-nadir  angle, 
y.  Then,  a  simple  physical  model  for  the  stationary  clutter  DC  (CDC)  can  be  written  in  the  general 
form: 

Mi)  =  frefii)  +  (<Min(y)  —  typCos(y))ffor  (1) 

where  fc  is  the  CDC  and  frej  and  fjor  arc  deterministic  parameters.  These  deterministic  values  can 
be  computed  as  a  function  of  the  off-nadir  angle  and  satellite  state  vectors  (SV’s),  which  include 
position  (rx,  ry  and  r~)  and  velocity  (vx,  vy  and  vz).  However,  SV,  as  well  as  attitude,  may  be 
expressed  in  any  one  of  the  several  widely-used  reference  systems.  Conventionally,  SV’s  arc  given 
either  in  the  Earth  Centred  Inertial  (ECI)  system  or  in  the  Earth  Centred  Earth  Fixed  (ECEF)  system. 
Attitude  is  measured  in  a  local  flight  system  that  can  also  be  defined  in  various  ways.  The  most 
common  four  reference  frames  are  considered  here.  Each  of  them  is  a  right-handed  system  with  axes 
pointing  in  the  forward,  right  and  downwards  directions,  relative  to  the  satellite’s  trajectory.  The 
downwards  direction  may  be  defined  in  the  geodetic  sense  (orthogonal  to  the  surface  of  the  Global 
Ellipsoid)  or  geocentric  sense  (pointing  to  the  Earth  center),  while  the  forward  direction  can  be 
defined  with  respect  to  the  satellite’s  motion  in  the  ECI  or  ECEF  system.  The  most  suitable  choice 
of  the  flight  reference  system  depends  on  the  characteristics  of  the  SAR  and  attitude  sensors.  The 
attitude  of  RADARSAT-1,  which  is  not  yaw  steered,  is  commonly  specified  in  the  geodetic,  ECI- 
related  local  system.  ENVISAT  is  yaw  steered  and  must  maintain  a  geodetically  horizontal  position, 
therefore  it  is  suitable  to  use  the  geodetic  local  system  related  to  the  ECEF  system.  RADARSAT-2 
is  also  yaw-steered,  but  its  reference  local  system  is  geocentric. 

The  following  set  of  calculations  accounts  for  all  of  these  varieties: 

0  for  geocentric 

{R2e/R-l  -  1  )(rz  —  H sin  9/a()  for  geodetic 

and 

0  for  ECEF 
coe  for  ECI 


(2) 

(3) 
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where  Re  and  Rp  arc  the  Earth  Global  Ellipsoid  semi-major  and  semi-minor  axes,  H  is  the  satel¬ 
lite  altitude,  8/at  is  the  geodetic  latitude  and  (Oe  is  the  Earth  rotation  speed  (assumed  constant 
7.2921 158553e-5  rad/s).  The  expression  for  8  is  valid  for  any  of  the  Global  Ellipsoid  models,  but 
we  arc  using  the  WGS84  model.  Using  these  values  and  the  corresponding  SV  (in  ECI  or  ECEF), 
further  calculations  arc  straightforward: 


r(8) 

=  \JrX  +  rj  +  (rz  +  §)2 

(4) 

Vyer(S) 

rxvx  +  ryVy  +  (rz  +  6)vz 
r(8) 

(5) 

Vhor(&) 

=  \Jvl  +  vj  +  vj-vler(  8) 

(6) 

Vfor{  8) 

rxvy-ryvx 

—  Vhorify  ® 

Vhor\P) 

(V) 

V  down($) 

=  -VVer(5) 

(8) 

^ right  (8) 

r(5)vz-(rz  +  5)vver(8) 

Vhor{ 8) 

(9) 

(10) 

For  a  SAR  operating  at  wavelength  X  (the  carrier  wavelength),  the  reference  values  in  (1)  arc: 


frefiy)  —  ^  (u/oH’n(8)  cosy-l-  vright  (8)  siny) 

ffor  =  —  j^Vyor(8) 


(11) 

(12) 


At  the  level  of  first  order  approximation,  CDC  depends  linearly  on  yaw  and  pitch  and  can  be  calcu¬ 
lated  easily  for  a  given  off-nadir  angle,  i.e.,  for  all  slant  range  bins. 


If  we  choose  the  geocentric,  ECEF  related  reference  frame  for  attitude,  RADARSAT-2’s  yaw¬ 
steering  is  included  in  /re/,  while  the  attitude  <j)v  and  in  (1)  represents  the  residual  offset  between 
the  actual  pointing  and  the  reference  (nominal)  pointing  based  on  the  yaw-steering  law.  For  this 
reason,  the  offset  can  be  considered  constant  over  a  period  of  slow  time.  In  particular,  if  (f>v  and  <\)f, 
arc  estimated  from  the  SAR  signal  at  some  place  along  the  orbit  that  has  suitable  clutter  properties, 
they  may  also  be  used  at  a  near-by  place  where  the  stationary  CDC  cannot  be  extracted  from  the 
SAR  signal  directly.  This  approach  is  used  to  establish  stationary  CDC  values  over  the  patches  of 
water  using  the  measurements  taken  over  the  land  within  the  same  SAR  scene. 


Unlike  RADARS  AT- 1 ,  RADARS  AT-2  auxiliary  data  provide  an  on-board  measurement  of  the  at¬ 
titude  angles  with  reasonable  accuracy.  According  to  RADARS  AT-2  documentation  [7],  the  auxil¬ 
iary  attitude  values  arc  in  the  geocentric,  ECI  related  reference  frame.  They  arc  variable  over  slow 
time  because  yaw  steering  is  time  variant  in  this  reference  frame.  The  auxiliary  attitude  values  arc 
recorded  (8  per  second)  and  can  be  used  in  the  CDC  calculations  if  direct  CDC  estimation  is  not 
possible.  This  is  an  alternative  approach  which  may  be  used  when  the  scene  is  entirely  over  ocean. 
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2.2  Beamforming  of  MODEX  data 


MODEX-l  mode  is  an  experimental  dual-channel  mode  designed  for  MODEX,  where  the  pulses 
arc  transmitted  by  the  full  antenna  and  received  in  parallel  by  the  two  wings  of  the  antenna,  usually 
referred  to  as  fore  and  aft.  This  configuration  can  be  used  for: 

•  Displaced  Phase  Center  Antenna  (DPCA)  technique  for  moving  target  detection; 

•  Along  track  interferometry  (ATI)  technique  for  estimation  of  target  radial  speed; 

•  Direction  of  Arrival  (DOA)  technique  for  estimating  the  azimuth  position  of  the  target. 

Apart  from  these  typical  dual-channel  applications,  the  fore  and  aft  signals  can  also  be  coherently 
combined  to  produce  the  equivalent  of  the  full  antenna  receive  signal.  Beamforming  is  a  linear 
coherent  operation  that  can  steer  the  beam  to  attain  a  maximum  (or  to  create  a  notch)  in  the  desired 
direction  (or  directions  when  the  number  of  channels  is  more  than  two).  In  our  case,  it  is  desirable 
to  restore  the  beam  pattern  of  the  original,  full  length  antenna  by  a  linear  combination: 

S  =  [SforeSaft}G  (13) 

where  G  is  the  so-called  steering  vector.  Assuming  that  beams  of  both  wings  point  in  the  same 
direction  and  assuming  that  the  beam  patterns  arc  identical  (in  magnitude  and  phase),  this  would 
be  achieved  by  simply  adding  the  fore  and  aft  signals,  i.e.,  GH  =  [11].  However,  the  fore  and 
aft  channels  have  a  slight  gain  imbalance  and  a  phase  offset,  which  needs  to  be  corrected  before 
summing  them  together.  Yet  another  problem  is  the  slight  pointing  difference  for  the  two  antenna 
beams,  which  is  evident  from  Doppler  centroid  estimates  performed  on  each  channel  separately. 
This  problem  could  be  mitigated  by  a  range-dependent  G,  but  the  range-dependent  effect  is  judged 
to  be  small  and  is  ignored  in  this  work.  It  is  expected  that  a  better  alignment  of  the  fore  and  aft 
beams  will  be  achieved  in  the  near  future  by  implementing  a  suitable  weighting  across  the  antenna 
columns.  This  correction  has  already  been  up-loaded  and  tested  for  some  modes,  but  has  not  been 
applied  to  the  beams  used  in  this  analysis.  Therefore,  there  is  a  small  loss  of  coherence  in  our  data 
sets. 


The  beamforming  operation  can  be  done  after  the  front-end  processing  and  subsequent  linear  oper¬ 
ations  in  the  SAR  processor.  In  this  case,  it  was  performed  on  the  PHP  after  range  compression. 

For  the  Gibraltar  data  sets  considered  here,  the  form  of  the  steering  vector  is: 


G 


Sc  exp(j<|)c) 


(14) 


and  the  optimum  gain  achieved  is  slightly  over  5  dB  (compared  to  the  fore  channel  alone).  Under 
ideal  conditions,  it  would  be  6  dB.  This  gain  is  measured  by  comparing  the  estimated  relative  a°  for 
some  selected  land  area  over  Spain  and  Morocco.  Alternative  steering,  such  as  broadside  (which 
would  create  zero  Doppler  conditions),  has  also  been  considered,  but  this  provided  a  lower  gain. 
At  the  same  time,  the  noise  level  increased  by  3  dB  compared  to  a  single  channel  acquisition  since 
noise  is  summed  incoherently  and  gc  «  1.  Therefore,  the  signal-to-noise  ratio  (SNR)  is  worse  than  in 
the  regular  single-channel  acquisition,  but  the  signal-to-clutter  ratio  (SCR)  is  preserved.  The  beam 
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width  corresponding  to  S fore  and  S aft  is  reduced  in  S  to  the  value  expected  for  the  full  antenna  and 
the  coherent  integration  time  is  reduced  accordingly. 

The  power  imbalance  compensation,  gc  is  easy  to  measure.  The  phase  imbalance  compensation, 
(j)c  can  be  found  by  searching  for  the  largest  attainable  beamforming  gain.  However,  a  more  direct 
method  was  used.  Based  on  many  data  sets,  it  has  been  found  that  the  average  ATI  phase  offset, 
§cai,  and  the  CDC  satisfy: 

cl 

tycal  —  27C  fc  (])c  (15) 

Vi 

where  d  is  the  effective  baseline  and  i’v  is  the  satellite  speed.  Proportionality  between  the  average 
ATI  offset  and  CDC  is  expected  because  both  depend  on  the  antenna  squint.  The  constant  term,  (])c 
describes  an  electronic  offset  between  the  channels.  Based  on  fitting  over  different  scenes,  its  value 
is  found  to  be  0.331  for  HH  polarization  and  -2.319  for  VV  polarization. 
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3  Gibraltar  data  sets 


The  data  sets  used  in  this  report  were  acquired  over  Gibraltar  in  the  MODEX- 1  mode,  with  proper¬ 
ties  summarized  in  Table  1  and  swaths  shown  in  Fig.  1 .  The  third  pass  is  the  same  as  the  second, 
but  taken  with  a  different  polarization  48  days  later.  All  data  sets  have  the  same  pulse  duration  of 
41 .6qs  and  the  same  pulse  bandwidth  (BW)  of  30  MHz,  which  are  almost  identical  to  the  pulse  pa¬ 
rameters  of  the  RADARS  AT- 1  Fine  Beam  mode.  The  pulse  repetition  frequency  (PRF)  is  somewhat 
higher  to  match  the  higher  Doppler  bandwidth  of  the  MODEX  signal. 


Table  1 :  Data  set  parameters. 


data  set 

date 

yyyy-mm-dd 

pass 

off-nadir 

angle 

[°] 

pulse 

BW 

[MHz] 

pulse 

polarization 

PRF 

[Hz] 

Dl 

2009-01-19 

ascending 

44.30 

30 

HH 

1802.09 

D2 

2009-01-21 

descending 

37.10 

30 

HH 

1825.15 

D3 

2009-03-10 

descending 

37.10 

30 

VV 

1825.15 

Figure  1 :  RADARSAT-2  swaths  considered  overlaid  on  Google  Map  (©2009  Google  -  Map  data 
©2009  Europa  Technologies). 
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3.1  Gibraltar,  2009-01-19 


A  map  of  the  local  CDC  was  produced  by  running  CHASP  on  small  image  chips  from  throughout 
the  entire  data  set.  This  scene  is  about  75  km  long.  A  land  mask  was  used  to  select  image  chips 
belonging  to  land.  After  CDC  fitting  over  land,  the  attitude  angles  were  estimated  to  be  0.0204° 
yaw  and  0.0192°  pitch,  which  is  within  RADARS AT-2  specifications  [8].  In  Fig.  2,  it  is  clear  that 
the  land  chips  correspond  to  much  higher  G°,  although  the  data  are  not  absolutely  calibrated.  This 
figure  also  shows  the  uniformity  of  the  scene,  which  is  important  for  good  DC  estimation.  In  the 
maps  shown  in  Fig.  3,  there  is  a  clear  difference  between  the  land  and  sea  surface.  The  surface 
currents  have  a  positive  radial  speed,  i.e.,  away  from  the  SAR,  which  means  West  to  East  for  this 
ascending,  right-looking,  scene.  There  appears  to  be  an  anomaly  in  the  CDC  maps  close  to  the 
coastlines.  There  is  a  reasonable  explanation  for  this  artifact.  It  may  be  attributed  to  the  sidelobes 
coming  from  the  land  and  contributing  to  the  estimated  CDC  in  the  areas  of  low  ocean  back  scatter, 
which  is  a  consequence  of  the  FTFT  polarization,  high  incidence  angle  and  relatively  calm  sea  surface. 

As  a  consistency  check,  the  CDC  map  was  compared  to  the  dual  channel  calibration  map,  estimated 
by  adaptively  balancing  the  fore  and  aft  channels,  and  also  by  running  CF1ASP  on  small  data  chips. 
Dual  channel  balancing  is  based  on  the  requirement  that  the  ATI  phase  should  be  zero  for  the 
stationary  clutter.  This  is  achieved  by  estimating  and  removing  the  phase  offset,  which  may  be 
present  between  the  co-registered  complex  images  from  the  two  channels.  This  phase  offset,  (j)ra/, 
may  be  spatially  variable,  just  like  CDC.  It  has  been  found  that  the  two  maps  are  highly  correlated 
(proportional,  but  with  opposite  signs).  This  is  expected  since  both  CDC  and  the  channel  phase 
imbalance  depend  on  the  antenna  squint.  This  is  an  important  feature  for  our  analysis  because  the 
beamforming  approach  depends  on  it.  Moreover,  the  same  artifact  occurs  along  the  shore  for  both 
methods. 


Figure  2:  Land  mask  and  relative  G°  map  for  the  entire  Gibraltar  scene  from  2009-01-19. 
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(b)  Residual  CDC  map 


Figure  3:  Local  calibration  maps  for  the  entire  Gibraltar  scene  from  2009-01-19. 


3.2  Gibraltar,  2009-01-21 

This  scene  is  around  155  km  long  and  includes  more  water  than  land.  For  the  entire  scene,  the  CDC 
map  is  compared  to  the  (j)ca/  map  as  shown  in  Fig.  4.  Both  maps  were  produced  by  CF1ASP  within 
identical  image  chips.  As  expected,  these  maps  are  correlated.  Outliers  appear  in  both  maps  and  can 
be  associated  with  the  local  intensity  contrast  and  motion.  There  is  also  a  clear  difference  between 
the  land  and  sea  surface.  The  surface  currents  have  a  negative  radial  speed,  i.e.,  towards  the  SAR, 
which,  again,  means  West  to  East  for  this  descending,  right-looking,  scene.  Certain  areas  in  this 
scene  have  very  large  CDC  deviations  compared  to  the  calibrated  physical  model.  Surface  motion 
in  such  areas,  if  uncompensated,  may  contribute  up  to  0.7  m/s  to  the  radial  speed  estimates.  Other 
areas  in  the  same  scene  are  less  affected.  There  are  no  noise-dominated  regions  in  this  scene.  For 
this  long  scene,  processing  was  done  in  three  separate  sub-scenes  and  more  refined  CDC  estimates 
were  produced  for  each.  The  first  sub-scene  is  over  Spain  and  served  for  CDC  calibration.  The 
results  of  local  CDC  estimation  are  presented  in  the  Fig.  5  for  the  heamformed  signal  and  compared 
to  the  predicted  DC  values,  derived  for  the  RADARS  AT- 2  nominal  steering  law.  In  this  case,  smaller 
image  tiles  are  used  to  produce  a  finer  map  than  in  Fig.  4  and  median  filtering  was  applied  within 
each  tile  to  minimize  the  effects  of  local  discontinuities  in  the  backscatter  reflectivity.  After  fitting, 
the  attitude  angles  are  estimated  to  be  0.0238°  yaw  and  0.0255°  pitch,  which  is  within  RADARS AT- 
2  specifications  and,  in  fact,  very  close  to  the  estimates  from  the  2009-01-19  scene.  Estimates  for 
separate  channels  were  somewhat  different,  as  discussed  earlier  (0.0464°  fore  yaw,  0.0403°  fore 
pitch,  0.0040°  aft  yaw,  0.0221°  aft  pitch).  The  fitting  error  is  shown  in  Fig.  6.  The  applied  CDC 
estimation  algorithm  works  best  for  uniform  scenes.  The  relative  0°  map  is  included  for  reference. 

The  calibrated  model  was  then  applied  to  the  other  two  sub-scenes  with  ocean  and  ships.  As 
discussed  earlier,  this  approach  relies  on  the  assumption  that  the  attitude  remains  the  same  in  the 
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Figure  4:  Local  calibration  maps  for  the  entire  Gibraltar  scene  from  2009-01-21. 


adopted  reference  frame.  In  other  words,  the  change  due  to  yaw-steering  is  accounted  for,  but  the 
deviation  from  it  is  assumed  to  be  constant.  This  assumption  is  now  confirmed.  Fig.  7  shows  the 
difference  between  the  calibrated  model  (from  Spain)  and  the  measured  CDC  for  the  second  sub¬ 
scene,  which  is  over  the  Atlantic  ocean,  with  no  land.  The  third  sub-scene  is  over  the  Atlantic  with 
the  Morocco  coast  at  near  range.  The  land  mask  and  the  relative  G°  map  are  shown  in  Fig.  8  and 
the  residual  CDC  is  shown  in  Fig.  9.  The  reference  attitude  for  this  sub-scene  is  still  from  the  Spain 
sub-scene,  but  does  match  the  CDC  well  over  the  land  of  Morocco. 

3.3  Gibraltar,  2009-03-10 

This  is  a  repeat  pass,  but  with  a  different  polarization.  The  data  were  acquired  48  days  (two  repeat 
pass  cycles)  after  the  2009-01-21  data.  CDC  was  calibrated  in  the  same  way,  using  only  the  patch  of 
data  over  Spain.  The  calibration  maps  are  shown  in  Fig.  10  and  11.  From  these  maps,  the  estimated 
attitude  for  this  scene  is  0.0308°  yaw  and  0.0236°  pitch.  All  ships  appear  in  a  sub-scene  over  the 
Atlantic.  The  G°  and  CDC  maps  for  the  Atlantic  sub-scene  are  shown  in  Fig.  12.  These  maps  are 
less  uniform  than  the  corresponding  maps  from  2009-01-21. 

3.4  Validated  ships 

There  are  17  AlS-validated  moving  ships  in  the  2009-01-19  data  set,  as  listed  in  Table  2.  There 
are  9  AlS-validated  moving  ships  in  the  2009-01-21  data  set,  as  listed  in  Table  3.  There  are  9  AIS- 
validated  moving  ships  in  the  2009-03-10  data  set,  as  listed  in  Table  4.  Their  length,  L,  speed,  V 
and  heading,  0  were  extracted  from  the  AIS  data.  The  radial  speed,  Vj,  can  be  derived  using  the 
imaging  geometry.  The  incidence  angle,  a ,-„c,  is  also  derived  from  the  geometry. 


DRDC  Ottawa  TM  2009-109 


9 


-0.132 

0.134 

-0.136 
-0.138 
-0.14 
-0.142  | 

9 

-0.144  ° 

-0.146 

-0.148 

-0.15 

-0.152 


(a)  CDC  map 


(b)  CDC  prediction  for  nominal  attitude 


Figure  5:  Local  calibration  maps  for  the  Spain  sub-scene  from  2009-01-21. 


Figure  6:  CDC  model  fitting  results  for  the  Spain  sub-scene  from  2009-01-21. 
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(a)  Calibrated  CDC  model 


(b)  Residual  CDC  map 


Figure  7:  Residual  CDC  for  the  Atlantic  sub-scene  from  2009-01-21. 


Figure  8:  Residual  CDC  for  the  sub-scene  off  Morocco  from  2009-01-21. 
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Figure  9:  Residual  CDC  for  the  sub-scene  off  Morocco  from  2009-01-21. 


Figure  10:  Land  mask  and  relative  a°  map  for  the  Spain  sub-scene  from  2009-03-10. 
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(b)  Residual  CDC  map 


Figure  11:  Local  calibration  maps  for  the  Spain  sub-scene  from  2009-03-10. 
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(a)  a°  map 


(b)  Residual  CDC  map 


Figure  12:  Local  calibration  maps  for  the  Atlantic  sub-scene  from  2009-03-10. 
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Table  2:  Ships  observed  on  2009-01-19. 


AIS 

derived 

ship 

L  [m] 

V  [m/s] 

0  [°T] 

Vt  [m/s] 

dine  [  ] 

D1S01 

86 

18.09 

45.9 

11.559 

51.7 

D1S02 

333 

6.06 

252.9 

-4.744 

52.3 

D1S03 

172 

6.43 

252.7 

-5.037 

52.4 

D1S04 

101 

5.14 

84.4 

4.050 

52.1 

D1S05 

99 

7.58 

76.1 

5.911 

51.5 

D1S06 

193 

6.38 

270.0 

-4.961 

51.9 

D1S07 

122 

6.93 

258.4 

-5.452 

52.0 

D1S08 

115 

9.44 

67.3 

7.203 

51.9 

D1S09 

87 

14.87 

246.8 

-11.332 

52.0 

D1S10 

181 

8.75 

271.9 

-6.708 

51.2 

dish 

115 

6.71 

89.0 

5.212 

51.6 

D1S12 

88 

3.64 

250.5 

-2.819 

52.2 

D1S13 

215 

7.48 

252.0 

-5.856 

52.6 

D1S14 

142 

7.56 

267.0 

-5.885 

51.9 

D1S15 

250 

6.45 

247.0 

-4.934 

52.1 

D1S16 

132 

4.94 

252.2 

-3.871 

52.5 

D1S17 

125 

5.30 

251.6 

-4.139 

52.5 

Table  3:  Ships  observed  on  2009-01-21. 


AIS 

derived 

ship 

L  [m] 

V  [m/s] 

0  [°T] 

Vt  [m/s] 

dine  [  ] 

D2S01 

147 

7.55 

212.6 

1.976 

42.6 

D2S02 

196 

9.92 

265.9 

6.489 

42.4 

D2S03 

123 

5.66 

229.2 

2.389 

41.9 

D2S04 

225 

5.33 

272.1 

3.566 

42.5 

D2S05 

69 

5.15 

89.5 

-3.382 

41.9 

D2S06 

200 

8.41 

225.7 

3.338 

42.6 

D2S07 

176 

7.67 

32.7 

-2.024 

42.8 

D2S08 

108 

4.67 

282.0 

3.148 

42.8 

D2S09 

157 

7.13 

71.2 

-4.148 

41.6 
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Table  4:  Ships  observed  on  2009-03-10. 


AIS 

derived 

ship 

L  [m] 

V  [m/s] 

0  [°T] 

Vt  [m/s] 

C^inc  [  ] 

D3S01 

128 

7.77 

271.1 

5.068 

41.3 

D3S02 

274 

7.60 

89.3 

-5.078 

42.8 

D3S03 

117 

4.97 

89.4 

-3.244 

41.6 

D3S04 

115 

6.10 

271.3 

4.004 

41.6 

D3S05 

205 

6.51 

52.7 

-2.973 

42.2 

D3S06 

190 

7.88 

93.4 

-5.315 

42.7 

D3S07 

225 

6.99 

271.7 

4.589 

41.5 

D3S08 

179 

6.43 

90.8 

-4.266 

42.2 

D3S09 

132 

6.61 

90.5 

-4.290 

41.2 
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4  Velocity  estimation  results 


Velocity  estimation  includes  radial  speed  estimation  and  heading  estimation.  Radial  speed  estima¬ 
tion  is  based  on  the  Doppler  offset.  Heading  estimation  is  based  on  estimation  of  the  aspect  angle 
of  the  ship  in  the  image  domain. 


4.1  Radial  speed  estimation  results 

Radial  speed  estimation  is  done  on  the  beamformed  PHP  using  the  same  Doppler  offset  algorithm 
and  software  as  was  used  in  the  case  of  RADARSAT-1  data  [2].  The  only  difference  is  that  local 
CDC  estimation  was  disabled  while  processing  ships.  A  new  utility  was  developed  to  provide  CDC 
values  according  to  the  calibrated  model,  which  were  fed  into  CHASP  as  external  settings.  With 
this  approach,  no  additional  corrections  are  needed. 


The  effect  of  beamforming  is  seen  in  Fig.  13  where  power  profiles  of  the  fore,  aft  and  beamformed 
signal  are  compared  at  the  output  of  the  adaptive  filter  used  for  SCR  and  SNR  enhancement  [3]. 
This  figure  illustrates  that  the  ship  signature  is  extracted  well  from  the  background  when  adaptive 
filtering  is  applied  to  each  of  the  channels  separately  and,  also,  when  applied  to  the  beamformed 
signal  in  slow  time.  The  sharper  slope  at  the  leading  and  trailing  end  of  the  ship’s  power  profile, 
after  beamforming,  makes  it  easier  to  estimate  the  acquisition  time  (beam  centre  time)  at  which  the 
Doppler  offset  must  be  estimated.  Results  for  2009-01-19  and  2009-01-21  are  summarized  in  the 
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Figure  13:  Power  envelope  extracted  by  the  Doppler  offset  algorithm. 


scatter  plot  in  Fig.  14.  The  over- all  standard  deviation  of  the  radial  speed  estimates  is  0.85  m/s, 
which  is  almost  the  same  as  the  accuracy  of  the  RADARSAT-1  estimates  obtained  from  the  Fine 
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Beam  scenes  over  the  Strait  of  Dover.  For  comparison,  the  standard  deviation  of  the  radial  speed 
estimates  based  on  the  dual-channel  ATI  method  is  0.3  m/s.  The  ATI  method  achieves  higher 
accuracy  at  the  expense  of  additional  receiver  hardware  and,  typically,  a  reduced  swath.  Both 
methods  face  the  same  problem  in  the  presence  of  surface  motion  since  they  both  require  calibration. 
Both  methods  gave  better  results  for  the  2009-01-19  scene  than  the  2009-01-21  scene.  In  the  latter 
scene,  the  incidence  angle  was  smaller  and  SCR  was  lower.  It  is  also  possible  that  the  weather 
conditions  contributed  to  the  difference  noted  between  these  two  scenes.  Radial  speed  estimation 


Figure  14:  Radial  speed  scatter  plot  for  the  validated  ships  in  HH  scenes. 

was  more  difficult  in  the  VV  scene  from  2009-03-10.  In  this  case,  the  adaptive  algorithm  for 
Doppler  tracking  was  faced  with  a  lower  SCR,  which  is  evident  in  Fig.  15.  As  a  result,  there 
was  a  larger  estimation  error.  As  demonstrated  in  Fig.  16,  the  speed  estimates  are  biased  towards 
zero.  An  explanation  for  this  bias  lies  in  the  nonstationary  character  of  the  clutter  upon  deramping 
so  that  tracking  initiates  under  different  clutter  conditions  depending  on  the  Doppler  offset  of  the 
tracked  ship  and  depending  on  the  starting  position  in  slow  time.  Due  to  this  effect,  there  is  an 
unequal  latency  of  the  forward  and  backward  tracking,  as  seen  in  Fig.  15.  For  a  ship  with  negative 
radial  speed  (positive  Doppler  offset),  forward  tracking  commences  slower  than  backward  tracking 
(case  (a)  in  Fig.  15).  For  a  ship  with  positive  radial  speed  (negative  Doppler  offset),  backward 
tracking  commences  slower  than  forward  tracking  (case  (b)  in  Fig.  15).  Being  non-symmetrical, 
tracking  delays  do  not  fully  compensate  each  other,  giving  rise  to  a  bias  in  the  estimate  of  the 
signature  acquisition  time,  which  affects  the  estimation  of  the  Doppler  offset.  The  same  effect  is 
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also  noticeable  in  Fig.  13  for  HH  data,  but  to  a  smaller  degree.  In  principle,  the  algorithm  could  be 
improved.  However,  possible  improvements  are  not  considered  in  this  report.  Standard  deviation 
of  the  estimation  error  was  roughly  50%  larger  than  for  the  2009-01-21  HH  scene  (acquired  in  the 
same  mode  over  the  same  region).  For  the  dual-channel  ATI  method,  standard  deviation  of  the  VV 
estimates  was  approximately  doubled  when  compared  to  the  standard  deviation  of  the  HH  estimates 
from  2009-01-21. 


azimuth  samples 


AN  F  forward 
—  ANF  backward 
-  ALE  forward 


azimuth  samples 


(a)  D3S08 


(b)  D3S01 


Figure  1 5 :  Power  envelopes  extracted  by  the  Doppler  offset  algorithm  in  the  VV  scene. 


4.2  Heading  estimation  results 

Relative  heading,  also  referred  to  as  the  aspect  angle,  can  be  estimated  in  the  image  domain. 

The  estimation  algorithm  based  on  region  growing  that  was  used  for  RADARS  AT- 1  data  [2]  was 
applied  in  this  case.  Compared  to  the  RADARSAT- 1  data  that  were  collected  in  the  Strait  of  Dover, 
there  were  more  estimation  problems  in  the  RADARS  AT- 2  Gibraltar  data  sets.  Different  problems 
occurred  for  different  scenes.  Therefore,  the  results  from  the  three  scenes  are  presented  separately. 

For  example,  the  results  for  2009-0 1-21  are  summarized  in  the  scatter  plot  in  Fig.  1 7  and  in  the  Table 
5.  Out  of  nine  ships  in  this  scene,  there  are  two  obvious  outliers,  D2S01  and  D2S07  and  yet  another 
ship,  D2S08,  has  a  large  error.  This  problem  was  explored  further  using  two  other  approaches  to 
estimate  the  orientation  of  the  ships.  Several  interesting  cases  were  studied  to  understand  better  the 
origin  of  the  problems. 

The  revised  algorithm  for  estimating  the  aspect  angle  starts  from  the  list  of  detected  ship  pixels. 
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Figure  16:  Radial  speed  scatter  plot  for  the  validated  ships  in  the  VV  scene. 


Assuming  that  (xk,yk),  k  £  (1,2, . . .  ,n),  represents  the  image  coordinates  of  a  detected  pixel  scaled 
to  the  units  of  length  on  the  ground  (meters),  the  matrix  of  the  detection  coordinates,  D  can  be 
formed  as  follows: 


D 


xi  y i 

*2  >’2 

Xn  yn 


(16) 


and  translated  to  the  centre  of  gravity: 


Df 


n 


\ 


i 

i 


(17) 


Ships  are  usually  elongated  in  one  direction,  so  that  the  Hermitian  matrix 

R  =  D^DCGM2  (18) 

tends  to  have  unequal  eigenvalues,  'kmax  and  /(„„■„ .  In  this  case,  the  major  eigenvector  \max  (cor¬ 
responding  to  Xmax)  indicates  the  dominant  orientation,  while  the  minor  eigenvector  Vm,„  is  its 
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Figure  17:  Aspect  angle  scatter  plot  for  the  validated  ships  on  2009-01-21. 


orthogonal  complement.  In  other  words,  the  rotation 

Q  =  \y  max  V/nin]  (19) 

reveals  the  major  and  minor  axes  of  the  vessel.  To  improve  accuracy,  points  in  D  are  weighted  by 
then-  corresponding  magnitudes  before  performing  the  eigenvalue  (EV)  decomposition.  An  illustra¬ 
tion  of  this  simple  concept  is  shown  in  the  example  of  two  ships  from  the  2009-01-21  data  set  in 
Fig.  18.  In  both  cases,  the  EV-based  orientation  is  visually  acceptable,  but  is  not  collinear  with  the 
AlS-derived  orientation.  In  most  cases,  the  EV-based  estimates  agree  well  with  the  estimates  from 
the  region  growing  algorithm,  as  presented  in  Table  5. 

The  discrepancy  between  our  estimates  and  the  AIS  data  may  be  explained  by  the  effects  described 
and  analysed  by  Ouchi  et  al.  [9].  In  the  wave  regime,  ships  exhibit  rotation  so  that  their  linear  speed 
varies  along  the  ship  (from  the  bow  to  the  stern)  and,  possibly,  with  height.  Thus,  rotational  motion 
can  cause  a  different  azimuth  shift  for  different  portions  of  the  ship,  which  produces  a  distortion  in 
the  image. 

To  further  investigate  this  assumption,  the  dual-channel  aspect  of  the  SAR  data  was  used.  The  two 
wings  of  the  antenna  form  a  two  element  array  that  can  be  used  to  reveal  the  direction  of  arrival. 
A  very  simple  DOA  algorithm  [10]  was  applied  to  all  pixels  to  estimate  their  position  in  azimuth. 
This  algorithm  is  not  sensitive  to  the  radial  speed  and  should  not  be  affected  by  ship  motion.  For 
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Table  5:  Aspect  angles  estimated  on  2009-01-21. 


algorithm 

region 

EV 

EV 

ship 

AIS 

growing 

detected 

located 

[°] 

[°] 

[°] 

[°] 

D2S01 

-67.3 

59.0 

61.9 

-75.4 

D2S02 

-13.9 

-17.5 

-5.3 

-18.0 

D2S03 

-50.7 

-74.3 

-75.8 

-66.5 

D2S04 

-7.6 

-35.8 

-29.0 

-5.6 

D2S05 

169.6 

154.9 

163.6 

103.8 

D2S06 

-54.1 

-71.8 

-74.3 

-77.8 

D2S07 

112.9 

-97.5 

-90.6 

94.2 

D2S08 

7.3 

-63.4 

-51.0 

-87.6 

D2S09 

151.3 

134.1 

135.3 

121.4 

-50 
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Figure  18:  Illustration  of  detected  ship  coordinates  with  the  major  EV  orientation  compared  to  the 
AlS-derived  orientation. 


the  detected  ship  pixels,  this  estimated  position  was  extracted  into  the  matrix  P,  similar  to  D.  The 
range  coordinates  in  P  and  D  are  the  same,  but  the  azimuth  coordinates  are  not.  In  P,  they  are 
estimated  via  the  DOA  algorithm,  while  in  D  they  are  the  apparent  locations  as  detected  in  the 
image.  The  same  EV  algorithm  was  applied  to  the  new  set  of  coordinates,  P,  as  was  applied  to 
D.  The  difference  is  illustrated  in  Fig.  19  for  two  ships.  In  both  cases,  the  aspect  angle  estimate  of 
the  located  ship  was  better  than  the  aspect  angle  estimate  of  the  detected  ship,  as  shown  in  Table  5. 
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However,  the  estimated  position  of  the  ship  is  not  very  crisp  due  to  DOA  estimation  errors,  which 
depend  on  the  SCR.  Therefore,  the  two-channel  DOA  method,  when  applied  to  each  pixel  in  the 
image  domain,  does  not  have  a  clear  advantage  over  the  single  channel,  image-based  method. 


Figure  19:  Illustration  of  detected  ship  coordinates  with  the  major  EV  orientation  compared  to  the 
estimated  ship  location  with  its  major  EV  orientation  and  AlS-derived  orientation. 


Finally,  the  images  of  the  two  outliers  are  shown  in  Fig.  20  and  overlaid  with  the  dual-channel 
processing  results.  The  apparent  orientation  of  D2S01  is  quite  different  from  the  DOA  result,  which 
is  in  agreement  with  the  AIS  data.  D2S07  is  a  nominally  along-track  traveling  ship  that  has  a  very 
elongated  signature  in  the  SAR  image.  It  is  interesting  to  note  that  the  DPCA  technique  detected 
two  separate  areas  of  the  ship  with  a  gap  in  between  (i.e.,  no  radial  motion)  and  the  ATI  technique 
found  both  positive  and  negative  radial  speeds.  The  DOA  algorithm  correctly  placed  the  detections 
within  the  gap  between  the  two  DPCA-detected  sections. 

We  may  conclude  that  the  wave  field  is  responsible  for  most  of  the  estimation  problems  in  the 
2009-01-21  data  set,  but  this  is  not  the  only  source  of  error. 

There  is  less  evidence  of  wave-induced  motion  in  the  2009-01-19  data  set.  In  this  image,  all  ships 
were  in  the  Strait  of  Gibraltar,  while  the  ships  from  2009-01-21  were  mostly  in  the  Atlantic  ocean. 
Almost  all  ships  in  the  2009-01-19  scene  were  range  traveling.  The  accuracy  of  the  aspect  angle 
estimation  for  these  ships  is  summarized  in  Fig.  21  for  the  EV  algorithm.  In  this  case,  there  are  also 
three  outliers,  all  passenger  crafts.  Two  of  them  are  very  fast  small  catamarans,  D1S01  and  D1S09, 
traveling  in  opposite  directions,  and  the  third  is  the  passenger  ship  D1S07.  Their  signatures  do  not 
have  the  expected  shape  in  the  image  domain. 

Several  estimation  problems,  caused  by  imaging  artefacts,  are  illustrated  below.  In  most  cases,  there 
is  a  smear  in  the  azimuth  direction.  Cargo  ship  D2S08  from  2009-01-21,  passenger  ship  D1S07 
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Figure  20:  Illustration  of  ship  imaging  in  a  wave  field  with  highlighted  pixels  detected  by  the  DPCA 
algorithm  and  overlaid  DOA-estimated  positions. 


from  2009-01-19,  and  small  passenger  catamarans  D1S01  and  D1S09  from  2009-01-19,  are  such 
examples.  It  is  difficult  to  speculate  on  the  source  of  these  effects.  For  example,  the  power  profile 
of  D2S08  demonstrates  a  sudden  strong  return  (Fig.  22)  that  could  have  been  caused  by  a  rotating 
reflective  surface  on  the  vessel.  For  D1S07,  a  similar  power  peak  coincides  with  a  brief  disruption 
in  the  Doppler  history  (Fig.  23),  which  supports  the  same  assumption.  Such  unusual  flashes  within 
the  vessel  signature  would  cause  a  smear  in  the  image  domain.  Images  of  D2S08  and  D1S07  are 
shown  in  Fig.  22  and  23.  Images  of  the  two  fast  crafts  are  shown  in  Fig.  24.  All  algorithms  for 
aspect  angle  estimation  failed  for  these  ships.  Although  small,  the  cargo  vessel  D1S12  (Fig.  25  (a)), 
does  not  have  the  same  problem  as  the  two  fast  crafts.  The  region  growing  algorithm  failed  in  the 
case  of  D1S06  from  2009-01-19  because  it  stopped  prematurely  due  to  a  discontinuity  in  the  ship’s 
signature  (Fig.  25  (b)).  The  EV  method  worked  well  in  this  case. 

In  the  2009-03-10  scene,  aspect  angle  estimation  was  mostly  good,  as  shown  in  Table  6  and  Fig.  26. 
Both  the  region  growing  algorithm  and  the  EV  algorithm  performed  well. 
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Figure  21:  Aspect  angle  scatter  plot  for  the  validated  ships  on  2009-01-19. 
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(a)  Power  profile  for  D2S08 


(b)  D2S08 


Figure  22:  A  short  strong  reflection  from  D2S08  may  have  caused  the  smear  in  the  azimuth  direc¬ 
tion. 


(a)  Doppler  history  for  D1S07 


Figure  23:  A  brief  disruption  in  the  Doppler  history  of  D1S07  may  be  due  to  the  interference  from 
a  moving  reflector  on  the  ship,  which  produced  the  azimuth  smear  causing  problems  in  aspect  angle 
estimation. 
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Table  6:  Aspect  angles  estimated  on  2009-03-10. 


algorithm 

ship 

AIS 

region 

growing 

EV 

detected 

[°] 

[°] 

[°] 

D3S01 

oo 

00 

28.4 

6.1 

D3S02 

169.5 

161.9 

169.4 

D3S03 

169.5 

146.4 

128.9 

D3S04 

-8.6 

-21.8 

-10.3 

D3S05 

132.8 

132.9 

127.6 

D3S06 

173.7 

159.7 

166.9 

D3S07 

-8.2 

-19.4 

-12.2 

D3S08 

171.0 

170.0 

173.3 

D3S09 

170.5 

154.4 

164.3 

Figure  26:  Aspect  angle  scatter  plot  for  the  validated  ships  on  2009-03-10. 
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4.3  Experiments  with  SLC  products 

Azimuth  focusing  influences  ship  signatures  in  the  following  ways: 

1.  The  Doppler  spectrum  is  usually  not  fully  preserved,  especially  if  spectral  windowing  is 
applied.  Since  the  Doppler  spectrum  of  a  ship  with  non-zero  radial  speed  has  a  frequency 
offset  relative  to  the  stationary  world  Doppler  spectrum,  spectral  shaping  is  non-synmietric 
for  the  ship.  A  portion  of  the  ship’s  Doppler  spectrum  may  be  cut  off.  This  would  adversely 
affect  the  Doppler  offset  algorithm. 

2.  Range  cell  migration  correction  (RCMC)  that  is  appropriate  for  the  stationary  world  is  not 
necessarily  suitable  for  moving  targets.  In  fact,  the  only  difference  is  in  the  positioning  of 
the  RCMC  discontinuity  in  the  Doppler  frequency  domain.  If  the  process  bandwidth  is  wide 
(equal  or  close  to  PRF),  a  portion  of  the  ship’s  signature  may  receive  an  inadequate  RCMC. 

3.  CDC  calibration  may  be  biased  if  the  Doppler  spectrum  is  shaped  according  to  the  internal 
estimates  of  the  SAR  processor.  Even  if  very  accurate  estimation  algorithms  arc  used  in  the 
SAR  processor,  there  can  be  a  bias  due  to  the  surface  motion. 

However,  if  the  PRF  is  high  and  the  processed  Doppler  bandwidth  is  large  and  well  preserved 
(i.e.,  a  boxcar  window),  the  SLC  products  may  be  a  good  starting  point  for  analysing  slow  moving 
vehicles  such  as  ships.  This  may  even  be  a  more  convenient  input  since  RCMC  would  not  be 
needed,  contrary  to  the  PHP-based  analysis. 

A  simple  experiment  was  devised  to  test  these  ideas.  The  PHP  data  were  first  focused  using  the 
Range-Doppler  algorithm  and  then  were  unfocused.  The  applied  Range-Doppler  focusing  algorithm 
included  all  standard  steps  and  the  processed  bandwidth  was  varied  from  full  PRF  to  70%  of  the 
PRF.  For  spectral  shaping,  CHASP  uses  the  Kaiser-Bessel  window  function  with  the  configurable 
shape  factor  [3.  Two  values  were  used,  [3  =  0  (boxcar  window)  and  p  =  2.8.  Unfocusing  was 
implemented  as  a  ID  operation  in  the  Doppler  frequency  domain  so  that  the  RCMC  operation  was 
not  undone.  The  main  effect  of  the  applied  procedure  is  the  controlled  reduction  of  the  Doppler 
bandwidth  and  the  implementation  of  the  RCMC,  based  on  the  stationary-world  DC. 

Reduction  of  the  processed  Doppler  BW  is  very  noticeable  in  the  power  envelope  of  the  adaptive 
filter  used  in  the  Doppler  offset  algorithm.  Fig.  27  (a)  represents  the  most  drastic  case  where  only 
70%  of  the  PRF  is  processed  in  the  case  of  the  fastest  ship  from  the  2009-01-19  data  set.  For 
better  clarity  of  the  plot,  only  the  forward  tracking  results  arc  shown.  Fig.  27  (b)  shows  the  case 
of  a  moderately  fast  ship  with  processed  BW  equal  to  90%  of  the  PRF.  Only  backward  tracking 
results  arc  shown  in  this  plot.  In  comparison  with  the  power  profiles  from  the  PHP  data,  the  SLC 
signatures  of  these  two  ships  arc  truncated  in  slow  time.  Although  this  truncation  is  much  less  for 
D2S02,  it  introduces  a  significant  bias  in  the  estimate  of  the  acquisition  time,  which  causes  a  bias  in 
the  Doppler  offset  estimate.  In  all  cases,  there  was  a  degradation  in  the  accuracy  of  the  radial  speed 
estimates.  Two  typical  plots  arc  shown  in  Fig.  28  for  a  relatively  slow  ship  and  for  a  moderately 
fast  one.  As  expected,  high-BW  SLC  products  appeal-  to  be  acceptable,  especially  for  slower  ships. 
Radial  speed  estimation  using  SLC  products  was  biased  for  most  cases.  Therefore,  various  SLC 
products  were  compared  using  the  root-mean-square  error  (rmse)  of  the  estimates  as  a  better  metric 
than  the  standard  deviation.  Tables  7  and  8  present  the  rmse  of  the  radial  speed  estimates  when  the 
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azimuth  samples  azimuth  samples 


(a)  D1S01  with  BW/PRF=0.7 


(b)  D2S02  with  BW/PRF=0.9 


Figure  27 :  Examples  of  power  profiles  extracted  from  PHP  and  SLC  products  for  two  ships. 


Figure  28:  Examples  of  accuracy  degradation  with  reduced  Doppler  bandwidth. 


processed  Doppler  BW  is  reduced.  Comparing  the  two  HH  scenes,  the  results  are  worse  for  the 
2009-01-19  scene  because  there  are  two  fast  crafts  present.  For  the  VV  scene,  the  rmse  is  larger, 
but  for  the  SLC  products  with  a  high  BW  (over  90%  of  the  PRF)  and  [3  =  0,  the  rmse  is  comparable 
to  the  rmse  for  the  PHP.  If  processed  BW  is  less  than  shown  in  the  table,  then  the  algorithm  fails  for 
many  ships.  As  expected,  spectral  shaping  with  p  >  0  leads  to  increased  estimation  errors. 
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Table  7:  Radial  speed  rmse  for  all  ships  in  the  SLC  products  with  [3  =  0. 


PHP 

SLC  | 

[3  =  0 

BW/PRF 

1.00 

1.00 

0.95 

0.90 

0.85 

2009-01-19  rmse  [m/s] 

0.8434 

1.2023 

2.5266 

3.6726 

4.4881 

2009-01-21  rmse  [m/s] 

0.8968 

0.9543 

1.0037 

1.2227 

1.4976 

2009-03-10  rmse  [m/s] 

1.3058 

1.5816 

1.6132 

1.7525 

2.4524 

Table  8:  Radial  speed  rmse  for  all  ships  in  the  SLC  products  with  [3  =  2.8. 


PHP 

SLC  [3  =  2.8 

BW/PRF 

1.00 

1.00 

0.95 

0.90 

2009-01-19  rmse  [m/s] 

0.8434 

3.3097 

3.8827 

4.1856 

2009-01-21  rmse  [m/s] 

0.8968 

3.9757 

3.9922 

4.0031 

2009-03-10  rmse  [m/s] 

1.3058 

2.4433 

2.5644 

2.8767 
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5  Conclusion 


Based  upon  the  analysis  of  35  known  ships  within  three  RADARS AT-2  scenes  acquired  in  MODEX 
mode  and  converted  to  PHR  the  following  conclusions  and  observations  are  offered: 

•  Adaptive  algorithms  developed  for  RADARS  AT- 1  are  applicable  to  RADARS  AT-2,  at  least 
when  using  modes  with  similar  resolution; 

•  The  accuracy  of  radial  speed  estimates  is  similar  to  that  obtained  with  RADARSAT-1  at 
similar  resolution; 

•  For  HH  scenes,  the  standard  deviation  of  the  radial  speed  estimates  was  0.85  m/s,  compared 
to  a  standard  deviation  of  0.3  m/s  achieved  by  ATI  in  the  special  MODEX  mode; 

•  The  adaptive  algorithm  worked  for  the  VV  scene  without  any  modification,  albeit  with  a 
larger  root-mean-square  error  of  1 .3  m/s,  compared  to  a  standard  deviation  of  0.6  m/s  achieved 
by  the  AT  algorithm; 

•  PHP  is  the  most  suitable  input  for  radial  speed  estimation  by  adaptive  analysis,  but  a  cus¬ 
tomized  SLC  product  is  a  reasonable  alternative; 

•  If  SLC  is  used  as  the  input,  it  is  important  to  ensure  a  high  enough  PRF  and  to  process 
practically  the  full  PRF  bandwidth  with  a  box-car  shaped  spectral  window; 

•  For  the  first  time,  VV  and  HH  results  are  compared  for  ship  radial  speed  estimation;  HH 
polarization  provided  a  higher  accuracy  (in  the  absence  of  ship  wake  analysis); 

•  The  importance  of  a  calibrated  physical  model  is  affirmed  for  radial  speed  estimation; 

•  Based  on  the  limited  data  set  available  in  this  analysis,  the  DC  and  attitude  properties  of 
RADARSAT-2  seem  to  be  stable  and  predictable  from  the  auxiliary  data,  which  may  extend 
the  applicability  of  the  physical  model  to  ocean  areas  without  nearby  land  available  for  cali¬ 
bration; 

•  A  new  approach  to  ship  aspect  angle  estimation  based  on  eigen  analysis  was  used  in  the  image 
domain; 

•  The  effect  of  wave-induced  motion  on  the  apparent  ship  aspect  angle  is  clearly  demonstrated, 
but  cannot  be  resolved  in  a  single  aperture  SAR  image  (if  no  wakes  arc  visible); 

•  Certain  image  artefacts  can  confuse  the  present  algorithm  for  aspect  angle  estimation,  al¬ 
though  this  problem  can  be  corrected  by  visual  inspection,  which  suggests  that  the  algorithm 
could  be  further  improved. 

In  an  operational  setting,  it  is  expected  that  SAR-derived  ship  signatures  will  be  fused  with  AIS 
data.  Any  left  over  ship  signatures  become  ships  of  interest  since  they  arc  either  non-compliant  or 
non-cooperative  with  respect  to  AIS  deployment.  It  is  pointed  out  that  these  ships  arc  of  greatest 
interest,  but  have  the  least  known  about  them.  As  such,  additional  processing  of  the  SAR  data 
to  derive  additional  ship  information  would  be  beneficial.  For  example,  single-beam  ship  velocity 
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estimation  algorithms  such  as  those  described  in  this  document  would  be  useful  to  better  forecast 
the  future  location  of  the  ships  of  interest  if  deploying  other  surveillance  assets. 
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Annex  A:  Updated  CHASP  configuration 


The  CHASP  configuration  file  is  documented  in  [1],  which  was  written  in  2006.  Since  that  time, 
CHASP  has  evolved.  New  functionality  was  added,  which  is  reflected  in  the  configuration  file.  The 
current  list  of  controls  that  can  be  set  via  the  configuration  file  is  described  in  full  for  completeness. 
All  keywords  and  field  values  are  in  ASCII  form.  Field  values  represent  inputs  that  are  internally 
interpreted  as  one  of: 

•  character  strings  (e.g.,  file  names); 

•  numbers  (i.e.,  integer  or  floating  point); 

•  flags  or  logical  values  (e.g..  Yes,  No) 

Some  fields  are  required  (mandatory),  some  arc  optional  and  some  arc  conditional,  as  indicated  in 
the  following  list. 

input_data_f  ile  (character  string)  provides  the  name  of  the  input  file  or  files  with  range-com¬ 
pressed  data.  This  entry  is  mandatory;  CHASP  will  fail  without  it.  For  single  channel  data, 
the  full  file  name  (and  path,  if  different  from  the  current  directory)  must  be  provided.  It  can 
be  any  string  and  may  have  an  arbitrary  extension.  For  two  channel  or  four  channel  data, 
only  the  root  paid  of  the  file  name  is  required.  The  completion  of  the  input  file  names  is 
automatic,  according  to  CHASP  conventions.  For  quad-channel  data,  the  completions  are: 
hh .  rc,  hv .  rc,  vv .  r c  and  vh .  rc.  For  dual-channel  data,  the  completions  are:  f  ore .  r c  and 
aft .  rc. 

prf_file  (character  string)  optionally  provides  the  name  of  the  input  file  with  SVs.  If  present, 
this  file  must  be  in  Matlab-4  format  and  must  contain  the  SAR  positions  in  the  ECEF  system 
at  the  rate  of  the  SAR  PRF,  starting  from  the  first  signal  pulse  (not  necessarily  the  first  range 
compressed  pulse).  Georeferencing  is  not  possible  if  this  file  is  not  present.  This  file  may 
also  contain  the  “PRF/V”  values  at  the  same  sampling  rate. 

channels  (integer  1,  2  or  4)  tells  CHASP  how  many  channels  there  are  in  the  data  set.  The  default 
value  is  1.  For  more  than  one  channel,  the  naming  conventions  apply. 

f  irst_line  (integer)  tells  CHASP  at  which  input  range  line  to  staid  processing,  counting  from  the 
beginning  of  the  input  range  compressed  file  (which  may  differ  from  the  staid  pulse  of  the 
signal  data).  The  default  is  from  staid  of  file. 

last_line  (integer)  tells  CHASP  where  to  stop  processing,  but  in  most  cases  it  has  no  effect. 
Since  CHASP  processes  one  and  only  one  block  of  data,  it  is  actually  the  block  size  that 
determines  the  last  processed  line. 

first_sample  (integer)  tells  CHASP  at  which  range  sample  to  staid  processing,  counting  from 
the  first  sample  of  the  range-compressed  data  (which  may  differ  from  the  start  sample  of  the 
signal  data).  The  default  is  from  staid  of  the  range-compressed  lines. 

last_sample  (integer)  tells  CHASP  what  the  last  sample  is  for  processing.  The  default  is  to  ingest 
and  process  to  the  end  of  range  compressed  lines. 
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azimuth_focus  (flag  Yes  or  No)  controls  whether  to  perform  azimuth  compression  or  not.  By  de¬ 
fault,  azimuth  compression  is  not  applied. 

ZD_frame  (flag  Yes  or  No)  is  used  conditionally  to  choose  between  the  zero-Doppler  and  beam- 
centre  frame  if  azimuth  compression  is  enabled.  It  has  no  effect  if  azimuth  compression  is 
disabled.  By  default,  data  arc  compressed  to  the  zero-Doppler  image  frame.  For  No,  they  arc 
processed  to  the  beam-centre  image  frame. 

RangeDoppler  (flag  Yes  or  No)  is  used  conditionally  to  specify  which  focusing  algorithm  to  use  if 
azimuth  compression  is  enabled.  If  not  selected,  the  Range-Doppler  algorithm  is  not  used  by 
default.  Other  possibilities  include  the  SPECAN  algorithm  and  the  rudimentary  ID  “matched 
filter”  algorithm  (inherited  from  an  earlier  airborne  SAR  processor,  which  can  tolerate  only 
small  unambiguous  DC  offsets  and  has  no  capability  for  RCMC  in  frequency  domain). 

specan  (flag  Yes  or  No)  is  used  conditionally  to  specify  which  focusing  algorithm  to  use  if  azimuth 
compression  is  enabled.  If  not  selected,  the  SPECAN  algorithm  is  not  used  by  default.  Other 
possibilities  include  the  Range-Doppler  algorithm  and  the  rudimentary  ID  “matched  filter” 
algorithm. 

azimuth_align  (flag  Yes  or  No)  is  used  conditionally  to  specify  if  any  azimuth  shift  should  be 
performed  when  the  number  of  channels  is  more  than  one.  This  option  can  be  used  for 
azimuth  co-registration.  If  alignment  is  desired,  it  is  done  by  interpolating  the  data  of  certain 
channels.  For  quad-channel  data,  it  is  usually  the  VV  and  VH  channels  that  have  to  be 
interpolated  to  match  the  acquisition  times  of  the  HH  and  HV  data.  For  dual-channel  data, 
the  aft  channel  might  be  shifted  in  slow  time.  The  default  is  not  to  do  any  azimuth  shifting. 
If  enabled,  the  amount  of  shifting  is  defined  in  other  fields. 

azimuth_shift_hh  (floating  point  number) 

azimuth_shift_hv  (floating  point  number) 

azimuth_shift_vv  (floating  point  number) 

azimuth_shift_vh  (floating  point  number)  defines  the  amount  of  azimuth  shift  for  different  chan¬ 
nels  and  arc  applied  if  azimuth  (slow  time)  shifting  is  enabled.  They  arc  specified  in  terms 
of  the  pulse  spacing.  For  quad-channel  data,  each  value  is  applied  to  the  named  polarization. 
For  dual  channel  data,  the  hv  value  is  applied  to  the  aft  channel.  These  values  have  no  effect 
on  the  single  channel  data.  The  default  values  arc  0. 

compensate_radial  (flag  Yes  or  No)  tells  CHASP  if  any  range  shifting  is  required  for  certain 
channels.  This  option  can  be  used  for  range  co-registration,  i.e.,  if  it  is  necessary  to  align 
delays  for  different  channels  and  if  relative  delays  arc  not  corrected  in  pre-processing.  The 
default  is  not  to  do  any  range  shifting.  If  enabled,  the  amount  of  range  or  fast  time  shifting  is 
defined  in  other  fields. 

range_shift_hh  (floating  point  number) 

range_shift_hv  (floating  point  number) 

range_shift_vv  (floating  point  number) 


36 


DRDC  Ottawa  TM  2009-109 


range_shift_vh  (floating  point  number)  defines  the  amount  of  range  shift  for  different  channels 
and  arc  applied  if  range  (fast  time)  shifting  is  enabled.  They  are  specified  in  terms  of  the 
range  spacing.  For  quad-channel  data,  each  value  is  applied  to  the  named  polarization.  For 
dual  channel  data,  the  hh  value  is  applied  to  the  fore  channel  and  the  hv  value  is  applied  to 
the  aft  channel.  For  single  channel  data,  the  hh  value  is  applied.  The  default  values  arc  0. 

rcmc  (flag  Yes  or  No)  is  used  to  indicate  if  RCMC  in  Doppler  domain  should  be  done  or  not. 
RCMC  is  always  included  in  the  Range-Doppler  algorithm,  but  it  may  be  desired  even  if 
data  are  not  focused. 

doppler_looks  (integer  1  or  2)  defines  how  many  Doppler  looks  to  produce  in  the  case  that  az¬ 
imuth  focusing  is  enabled.  The  default  value  is  1 .  There  arc  other  fields  that  define  what  to 
do  with  looks  if  two  of  them  are  created. 

window_shape  (floating  point  number)  defines  the  shape  of  the  Kaiser-Bessel  window  if  filters 
arc  applied  in  either  range  or  azimuth.  The  default  value  is  0.0,  which  means  a  box-car 
window. 

full_bandwidth  (floating  point  number)  specifies  the  available  DBW  normalized  by  the  PRF  if 
azimuth  filtering  (focusing  or  shifting)  is  applied.  If  focusing  is  enabled,  this  field  must  be 
defined. 

look_bandwidth  (floating  point  number)  specifies  the  processed  DBW  per  look  normalized  by 
the  PRF  if  azimuth  filtering  (focusing  or  shifting)  is  applied.  If  focusing  is  enabled,  this  field 
must  be  defined  and  not  be  larger  than  the  full  bandwidth. 

look.of  f  set  (floating  point  number)  tells  CHASP  how  to  position  Doppler  looks  relative  to  the 
DC  and  is  expressed  in  fractions  of  the  PRF.  The  default  value  is  0.  For  non- zero  values,  the 
created  Doppler  look  or  looks  would  not  be  symmetric  about  the  DC. 

look.overlap  (floating  point  number)  specifies  the  amount  of  overlap  between  looks  as  a  frac¬ 
tion  of  the  PRF.  It  is  applicable  only  if  two  looks  arc  created. 

coherence  (flag  Yes  or  No)  tells  CHASP  if  coherence  between  azimuth  looks  should  be  com¬ 
puted.  This  option  makes  sense  only  under  certain  settings  that  arc  explained  in  detail  in 
[1]  (when  two  complex  azimuth  looks  are  being  created  and  azimuth  resizing  is  enabled). 
The  default  value  is  No. 

mag_coherence  (flag  Yes  or  No)  is  used  conditionally  to  choose  between  different  ways  to  com¬ 
pute  coherence.  The  default  value  is  No,  which  means  that  coherence  is  measured  on  complex 
looks,  while  Yes  tells  CHASP  to  measure  coherence  on  the  magnitude  of  the  azimuth  looks. 

norm.coherence  (flag  Yes  or  No)  is  used  conditionally  to  choose  between  different  ways  to  com¬ 
pute  coherence.  By  definition,  coherence  is  a  normalized  value,  but  this  is  not  the  default 
in  CHASP.  The  default  is  No,  which  produces  a  spatially  smoothed  interferogram  between 
azimuth  looks. 

multi_look  (flag  Yes  or  No)  controls  whether  processed  data  should  be  converted  to  magnitude, 
thus  losing  the  phase  information.  The  default  is  to  leave  data  as  complex  values.  This  flag 
controls  single  look  and  two  look  data,  focused  or  not. 
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spatial_size  (flag  Yes  or  No)  tells  CHASP  whether  resampling  in  azimuth  is  needed  to  produce 
square  pixels  of  ground  coverage.  There  is  no  provision  to  control  the  square  pixel  size.  The 
data  can  only  be  resampled  with  smoothing  in  azimuth  to  match  the  spacing  in  ground  range. 
The  same  operation  is  used  in  the  evaluation  of  inter-look  coherence. 

correlate_looks_len  (integer)  is  used  conditionally  to  specify  the  length  of  the  cross-correlation 
in  the  azimuth  direction  between  Doppler  looks.  It  is  used  only  if  two  looks  arc  created  and 
converted  to  magnitude.  Its  purpose  is  to  estimate  target  misregistration  in  the  two  looks.  The 
default  value  is  0  (not  applied). 

auto_fit  (flag  Yes  or  No)  is  used  to  enable  polynomial  fitting  of  the  Doppler  phase  modulation 
revealed  by  measuring  look-to-look  misregistration  (as  paid  of  the  “map-drift”  algorithm). 
The  default  value  is  No.  This  feature  can  be  used  only  iteratively,  i.e.,  CHASP  has  to  run 
multiple  times  with  positive  correlate_looks_len  and  it  can  perform  adaptive  fitting  in  the 
subsequent  run. 

range_look  (flag  No  or  Up  or  Low)  provides  the  option  of  creating  a  specified  range  look.  The 
default  option  is  No,  which  leaves  the  range  spectrum  unchanged.  The  other  two  options  arc 
used  to  extract  the  upper  or  the  lower  band  from  the  full  range  BW  (i.e.,  pulse  BW). 

range_bb  (flag  Yes  or  No)  is  used  conditionally  to  reposition  (in  frequency)  the  selected  range  look 
to  the  base  band.  This  is  a  step  in  range  multi-looking  and,  also,  a  step  in  some  algorithms 
for  PRF  ambiguity  resolution.  However,  CHASP  does  not  have  full  capabilities  for  such 
algorithms  (they  would  have  to  be  completed  externally). 

chirp_bw  (floating  point  number)  defines  the  pulse  BW  normalized  by  the  sampling  frequency 
and  serves  for  extraction  of  the  range  looks,  if  required.  Otherwise,  it  has  no  effect. 

doppler_grid_rg  (integer)  sets  the  number  of  grid  elements  in  the  range  direction  for  calculating 
the  DC  estimates.  DC  estimation  is  disabled  if  this  number  is  0,  which  is  the  default. 

geoid_over_ellipsoid  (floating  point  number)  is  used  to  specify  the  terrain  elevation  over  the 
reference  Global  Ellipsoid  (WGS84).  It  is  used  for  the  purpose  of  georeferencing,  if  applied. 
The  default  value  is  0.0. 

spread_loss  (integer)  is  the  degree  of  the  power  spread  loss  function  to  be  applied  in  range.  The 
default  is  0,  which  makes  no  radiometric  changes  to  the  data.  A  positive  number  can  be  used 
to  correct  possible  deficiencies  of  the  radiometric  normalization  applied  in  the  pre-processing. 

ref  erence_range  (floating  point  number)  is  used  as  a  reference  to  normalize  the  slant  range  be¬ 
fore  evaluating  the  spread  loss  radiometric  correction,  if  required. 

Pauli  (flag  Yes  or  No)  tells  CHASP  if  the  conversion  to  the  Pauli  polarimetric  basis  should  take 
place.  This  option  is  valid  only  for  quad-channel  data  on  condition  that  all  four  channels  arc 
processed  concurrently,  i.e.,  loaded  together  into  CHASP’s  memory  space.  The  default  value 

is  No. 

DPCA  (flag  Yes  or  No)  is  used  to  find  the  channel  difference,  sample  by  sample.  This  option  is  valid 
only  for  dual  channel  data  on  condition  that  both  channels  arc  processed  concurrently,  i.e., 
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loaded  together  into  CHASP’s  memory  space.  The  default  value  is  No.  This  flag  controls 
only  the  subtraction  function  and  no  other  step  of  the  standard  DPCA  technique,  which  is 
controlled  via  other  settings  (e.g.,  channel  co-registration  and  balancing).  Therefore,  this  flag 
can  be  useful  for  other  operations  such  as  beamforming. 

ATI  (flag  Yes  or  No)  is  used  to  find  the  channel  interferogram,  sample  by  sample.  This  option  is 
valid  only  for  dual  channel  data  on  condition  that  both  channels  arc  processed  concurrently, 
i.e.,  loaded  together  into  CHASP’s  memory  space.  The  default  value  is  No.  This  flag  controls 
only  the  complex  conjugate  multiplication  function  and  no  other  step  of  the  standard  ATI 
technique,  which  is  controlled  via  other  settings  (e.g.,  channel  co-registration  and  balancing). 
Therefore,  this  flag  can  be  useful  for  other  operations  such  as  the  DOA  algorithm. 

DOA  (flag  Yes  or  No)  is  used  to  enable  the  essential  step  of  the  implemented  DOA  algorithm  and 
that  is  the  multiplication  of  the  aft  channel  with  a  ramp  phase  function  in  slow  time.  This 
option  is  valid  only  for  dual  channel  data  on  condition  that  both  channels  are  processed  con¬ 
currently,  i.e.,  loaded  together  into  CHASP’s  memory  space  and  the  effective  ATI  baseline 
must  be  defined.  The  default  value  is  No.  This  flag  is  responsible  only  for  the  de-ramping  step 
of  the  DOA  algorithm  (ATI  flag  must  be  set  to  complete  the  algorithm).  If  applicable,  there 
arc  separate  optional  fields  to  specify  the  reference  point  for  de -ramping. 

ADI  (flag  Yes  or  No)  is  used  to  enable  the  essential  step  of  the  “azimuth  displacement  indicator”, 
an  algorithm  that  achieves  the  same  functionality  as  the  standard  ATI  algorithm.  This  option 
is  valid  only  for  dual  channel  data  on  condition  that  both  channels  arc  processed  concur¬ 
rently,  i.e.,  loaded  together  into  CHASP’s  memory  space  and  the  effective  ATI  baseline  must 
be  defined.  The  default  value  is  No.  This  flag  is  responsible  only  for  the  multiplication 
of  the  aft  channel  with  a  ramp  phase  function  before  and  after  azimuth  focusing  (ATI  and 
azimuth.focus  flags  must  be  set  to  complete  the  algorithm).  If  applicable,  there  arc  separate 
optional  fields  to  specify  the  reference  point  for  de -ramping. 

balance_dpca  (flag  Yes  or  No)  is  used  to  enable  channel  balancing.  This  option  is  valid  only  for 
dual  channel  data  on  condition  that  both  channels  arc  processed  concurrently,  i.e.,  loaded 
together  into  CHASP’s  memory  space.  The  default  value  is  No.  If  this  flag  is  set,  channel 
imbalance  is  evaluated  (including  the  ratio  of  the  average  magnitude  and  the  average  of  the 
phase  difference)  and  is  compensated  by  multiplication  of  the  aft  channel  with  a  complex 
scalar.  However,  the  estimated  imbalance  may  be  overwritten  by  external  values. 

balance_ratio  (floating  point  number)  sets  an  external  value  to  overwrite  CHASP’s  internal 
measurement  of  the  channel  magnitude  imbalance.  The  corresponding  compensation  is  ap¬ 
plied  to  dual  channel  data  if  this  value  is  positive  and  if  channel  balancing  is  enabled. 

balance_rad  (floating  point  number)  sets  an  external  value  to  overwrite  CHASP’s  internal  mea¬ 
surement  of  the  channel  phase  imbalance.  This  value  is  given  in  radians.  The  corresponding 
compensation  is  applied  to  dual  channel  data  if  balance-ratio  has  a  positive  value  and  if 
channel  balancing  is  enabled.  The  default  value  is  0.0.  An  external  value  can  be  used  if  any 
lineal-  combination  of  two  channels  is  desired,  e.g.,  for  beamforming. 

auto_focus  (flag  Yes  or  No)  is  used  to  enable  the  Doppler  frequency  tracking  algorithm.  Addi¬ 
tional  parameters  related  to  this  algorithm  can  be  set  through  other  fields.  The  default  value 
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is  No.  This  algorithm  cannot  run  together  with  azimuth  focusing.  If  tracking  is  enabled, 
focusing  is  disabled  automatically. 

alpha  (floating  point  number)  is  a  design  parameter  for  the  Doppler  tracking  algorithm  used  if 
auto_focus  is  enabled.  It  controls  the  initial  notch  width  of  the  adaptive  filter.  The  recom¬ 
mended  value  is  0.8,  which  is  also  the  default. 

track.range  (integer)  is  optionally  used  to  indicate  the  range  bin  at  which  the  target  of  interest 
is  located.  If  provided,  this  input  is  used  by  the  tracking  algorithm  and  for  extracting  various 
target  features  on  condition  that  such  algorithms  arc  enabled.  This  parameter  is  not  neces¬ 
sarily  used  in  a  strict  way  since  the  corresponding  CHASP  algorithms  may  search  for  the 
brightest  point  in  the  neighbourhood  of  the  given  range  bin.  If  not  provided,  these  algorithms 
will  pick  the  brightest  focused  point. 

track.azim  (integer)  is  optionally  used  to  indicate  the  azimuth  bin  at  which  the  target  of  interest  is 
located.  If  provided,  this  input  is  used  by  the  tracking  algorithm  and  for  extracting  various  tar¬ 
get  features  on  condition  that  such  algorithms  are  enabled.  This  parameter  is  not  necessarily 
used  in  a  strict  way  since  the  corresponding  CHASP  algorithms  may  search  for  the  brightest 
point  in  the  neighbourhood  of  the  given  azimuth  bin.  If  not  provided,  these  algorithms  will 
pick  the  brightest  focused  point. 

bunch.r  (integer)  is  optionally  used  to  indicate  the  number  of  range  bins  occupied  by  the  target  of 
interest.  If  provided,  this  input  is  used  by  the  tracking  algorithm  and  for  extracting  various 
target  features.  If  not  provided,  a  suitable  value  is  found  automatically  by  the  corresponding 
CHASP  algorithms. 

bunch_a  (integer)  is  optionally  used  to  indicate  the  number  of  azimuth  bins  occupied  by  the  target 
of  interest.  If  provided,  this  input  is  used  by  the  tracking  algorithm  and  for  extracting  various 
target  features.  If  not  provided,  a  suitable  value  is  found  automatically  by  the  corresponding 
CHASP  algorithms. 

track_strict  (flag  Yes  or  No)  is  used  to  enforce  the  use  of  track_range  and  track_azim  in  a 

strict  way,  i.e.,  without  any  adjustments  that  CHASP  algorithms  would  otherwise  apply.  The 
default  value  is  No. 

phase.order  (integer  0  to  4)  defines  the  polynomial  order  used  to  fit  the  phase  modulation  re¬ 
vealed  by  Doppler  tracking  or  by  iterative  use  of  look  correlation. 

auto  correct  (flag  Yes  or  No)  is  used  to  enable  adaptive  phase  compensation  in  slow  time  do¬ 
main  using  the  estimation  results  from  previous  CHASP  runs.  The  default  value  is  No.  If 
enabled,  this  feature  uses  either  a  polynomial  phase  compensation  (if  phase_order  is  posi¬ 
tive)  or  a  non-parametric  phase  compensation  extracted  previously  by  the  tracking  algorithm 
and  smoothed. 

smooth  (integer)  defines  the  size  of  the  weighted  moving  average  filter  applied  to  smoothen  the 
tracking  results  if  auto_correct  is  applied.  This  value  is  used  both  for  the  extracted  Doppler 
phase  variations  and  for  the  extracted  range  variations  of  the  tracked  target  over  slow  time. 
The  recommended  value  is  20,  which  is  also  the  default. 
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ACHIRP  (flag  Yes  or  No)  plays  a  similar  role  as  auto_correct,  but  applies  a  second  order  phase 
compensation  defined  by  the  along-track  speed,  which  is  set  in  a  separate  field. 

range_migration  (flag  Yes  or  No)  is  used  to  enable  RCMC  in  the  time  domain.  The  default  value 
is  No.  This  option  is  typically,  but  not  automatically,  selected  if  tracking  or  SPECAN  algo¬ 
rithms  are  used.  If  applicable,  there  arc  separate  optional  fields  to  specify  the  reference  point 
for  RCMC  in  the  time  domain. 

adaptive-migration  (flag  Yes  or  No)  is  used  to  enable  adaptive  RCMC  in  the  time  domain.  The 
default  value  is  No.  This  feature  can  be  used  only  iteratively,  i.e.,  CHASP  has  to  run  once 
with  auto_focus  enabled  and  can  apply  adaptive  RCMC  in  the  next  run. 

skew.az im  (integer)  is  used  conditionally  to  specify  the  reference  pulse  position  for  RCMC  in  time 
domain,  DOA  or  ADI  de-ramping.  The  default  value  is  the  centre  of  the  processing  block. 

locate  (flag  Yes  or  No)  is  used  to  request  calculation  of  the  target  radar  cross  section  (RCS)  and 
clutter  G°  and,  possibly,  other  target  attributes.  The  default  value  is  No.  If  data  are  converted 
to  magnitude  (multi -look  option  turned  on),  this  flag  enables  the  region  growing  algorithm, 
which  estimates  the  extent  and  orientation  of  the  target.  Other  attributes  include  the  slant 
range,  time,  ambiguity  distance  in  range  and  azimuth  and  an  estimate  of  the  ambiguity-to- 
target  ratio.  If  the  “prf  ’  file  is  provided,  incidence  angle  and  geographic  coordinates  of  the 
target  arc  calculated,  as  well.  By  default,  the  target  is  defined  as  the  brightest  spot  in  the  image 
with  its  immediate  surrounding  (as  determined  by  the  region  growing  algorithm).  However, 
if  track.range  and  track_azim,  bunch.r  and  bunch_a  entries  are  provided,  they  limit  the 
area  in  which  CHASP  searches  for  a  maximum,  while  track-strict  prevents  any  search 
and  enforces  the  external  definition  of  the  target. 

measure_contrast  (flag  Yes  or  No)  is  used  to  request  the  evaluation  of  contrast  in  the  area  around 
the  target.  The  default  value  is  No.  This  feature  can  be  used  only  if  locate  is  enabled. 
The  area  in  which  contrast  is  measured  is  automatically  defined  based  on  the  target  bounds 
as  estimated  by  the  region  growing  algorithm,  but  it  can  also  be  set  to  a  desired  value  via 
separate  optional  fields. 

window_range  (integer)  is  optionally  used  to  indicate  the  number  of  range  bins,  centred  at  the  tar¬ 
get  position,  where  contrast  must  be  measured  if  measure_contrast  is  enabled.  If  omitted, 
an  adaptive  window  is  used  based  on  the  region  growing  algorithm. 

window_azimuth  (integer)  is  optionally  used  to  indicate  the  number  of  azimuth  bins,  centred  at 
the  target  position,  where  contrast  must  be  measured  if  measure_contrast  is  enabled.  If 
omitted,  an  adaptive  window  is  used  based  on  the  region  growing  algorithm. 

profile  (flag  Yes  or  No)  is  used  to  request  the  evaluation  of  the  power  profile  along  azimuth  by 
averaging  in  range.  The  result  is  expressed  in  the  dB  scale.  The  number  of  averaged  range 
bins  depends  on  the  range  extent  of  the  target  (as  estimated  or  set  externally)  and  on  the  range 
displacement  of  the  ambiguities.  This  feature  can  be  used  for  plotting  and  visualizing  the 
ambiguities.  The  default  value  is  No. 

pgm.out  (flag  Yes  or  No)  is  used  to  create  a  PGM  file  as  CHASP  output.  The  default  value  is  No. 
This  option  is  valid  if  multi-look  is  on. 
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ppm.out  (flag  Yes,  No  or  Pauli)  is  used  to  create  a  PPM  tile  as  CHASP  output  where  channels 
arc  colour  coded.  The  default  value  is  No.  This  option  is  valid  only  for  quad-channel  data  if 
multi_look  is  on.  The  Pauli  option  is  allowed  only  if  all  channels  arc  processed  concur¬ 
rently,  i.e.,  loaded  together  into  the  CHASP’s  memory  space. 

pgm.floor  (floating  point  number)  is  used  conditionally  to  specify  the  cut-off  magnitude  in  dB 
(relative  to  the  maximum  intensity)  if  CHASP  is  doing  conversion  to  an  8-bit  output  format 
such  as  PGM  and  PPM.  This  feature  is  useful  for  creating  a  preview  of  the  processed  data. 

relative_dc_of  f  set  (floating  point  number)  sets  the  DC  value  normalized  by  the  PRF  that  is  to 
be  used  in  many  algorithms  such  as  focusing,  RCMC  or  azimuth  shifting.  The  default  value 
is  0.0. 

relative_dc_clutter  (floating  point  number)  sets  the  DC  value  normalized  by  the  PRF  that  is 
representative  of  the  background  clutter.  This  value  is  not  used  for  processing.  It  is  used  as 
a  reference  for  radial  speed  estimation  based  on  the  difference  between  the  estimated  target 
DC  and  this  clutter  value.  The  default  value  is  0.0. 

relative_v_of  f  set  (floating  point  number)  sets  the  target  along-track  speed  normalized  by  the 
SAR  speed.  It  can  be  used  to  modify  the  Doppler  rate  for  azimuth  focusing.  Alternatively,  it 
can  be  used  to  compensate  the  effect  of  the  target  along-track  speed  if  ACHIRP  is  used.  The 
default  value  is  0.0.  This  feature  is  often  used  in  iterative  CHASP  runs,  where  an  estimated 
value  from  one  run  becomes  an  input  parameter  for  the  next. 

prf  over  v  (floating  point  number)  is  used  to  optionally  enter  the  focusing  parameter  fp/vejj, 
which  is  the  ratio  of  the  PRF  and  the  effective  speed  of  the  SAR  platform.  This  parameter  is 
needed  if  azimuth_focus  is  enabled.  Another  possible  source  of  this  parameter  is  the  “prf” 
file,  but  can  be  overwritten  by  this  configuration  field.  If  neither  source  of  fp /veff  is  available, 
CHASP  will  fail. 

sin_inci_ang  (floating  point  number)  is  used  to  optionally  enter  the  sine  of  the  incidence  angle, 
which  is  needed  for  some  operations  (spatial_size  or  calculation  of  a°).  Another  possible 
source  of  this  parameter  is  the  “prf”  file,  but  this  is  an  option  to  provide  the  incidence  angle 
(externally  calculated)  even  if  the  “prf”  file  is  missing. 

wavelength  (floating  point  number)  is  the  only  way  to  define  the  SAR  wavelength  for  CHASP. 
It  is  required  for  many  calculations  and  most  functions  will  fail  if  it  is  not  provided. 

range_slant_spacing  (floating  point  number)  is  the  only  way  to  define  the  slant  range  sample 
size  for  CHASP.  It  is  required  for  many  calculations  and  most  functions  will  fail  if  it  is  not 
provided. 

azimuth_beamwidth  (floating  point  number)  provides  the  azimuth  beam  width  in  degrees,  which 
is  used  to  determine  the  length  of  the  tracking  window  when  auto_focus  is  enabled.  It  is 
required  for  the  tracking  algorithm. 

baselineATI  (floating  point  number)  provides  the  effective  baseline  length  that  is  used  if  D0A  or 
ADI  are  enabled.  It  is  not  needed  otherwise. 
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block.lines  (integer)  defines  the  size  of  the  data  block  in  azimuth.  Currently,  the  maximum 
allowed  value  is  16  K.  In  the  current  implementation,  this  value  does  not  have  to  be  a  power 
of  2.  The  default  value  is  1  K,  which  is  insufficient  if  CHASP  is  used  for  RADARSAT 
focusing  (CHASP  would  fail). 

processors  (integer)  defines  the  number  of  parallel  processes  that  will  be  used  to  execute  CHASP. 
The  default  value  is  1.  Almost  all  functions  can  be  executed  by  parallel  multiple  tasks.  Some 
adaptive  functions,  namely  tracking  and  region  growing,  arc  not  designed  to  be  parallelized. 

parts  (integer)  defines  the  number  of  parts  per  block  of  each  channel  when  CHASP  does  data 
partitioning.  The  default  value  is  1.  If  the  number  of  parts  and  the  number  of  processors  are 
both  greater  than  1,  different  processors  will  be  in  charge  of  different  parts. 

channels -together  (integer)  defines  the  number  of  channels  to  be  loaded  together  into  CHASP’s 
memory  space  and  processed  concurrently.  The  default  value  is  1.  For  single  channel  data, 
this  number  must  be  1.  For  dual  channel  data,  this  number  can  be  1  or  2,  but  certain  functions 
require  the  presence  of  both  channels  (e.g.,  DPCA,  ATI).  For  quad-channel  data,  this  number 
can  be  1,  2  or  4,  but  certain  functions  require  the  presence  of  all  channels  (Pauli). 

interp_size  (integer)  defines  the  length  of  the  “sine”  interpolation  kernel  used  for  resampling. 
In  the  current  implementation,  it  must  be  an  odd  number.  The  default  value  is  5,  but  the  usual 
value  is  1 1 . 

interp_f ractions  (integer)  defines  the  granularity  of  the  sub-pixel  shifts,  i.e.,  the  number  of 
different  “sine”  interpolation  kernels  used  for  resampling.  The  default  value  is  16,  but  the 
usual  value  is  64. 

change_input_endian  (flag  Yes  or  No)  is  used  to  request  byte  swapping  at  input.  The  default 
value  is  No.  This  option  can  be  useful  when  processing  data  created  with  a  different  byte 
order. 

change_output_endian  (flag  Yes  or  No)  is  used  to  request  byte  swapping  at  output.  The  default 
value  is  No.  This  option  can  be  useful  when  processing  data  that  will  be  used  on  some  platform 
with  a  different  byte  order. 
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Annex  B:  Data  processing  procedure 


The  Gibraltar  scenes  were  processed  systematically  using  a  slightly  modified  version  of  the  proce¬ 
dure  that  had  been  applied  earlier  to  the  RADARS  AT- 1  Dover  scenes.  The  difference  is  only  in  the 
pre-processing  and  data  preparation,  while  the  actual  CHASP  algorithms  arc  the  same. 

The  processing  procedure  is  as  follows: 

1.  Data  preparation,  which  includes  two  steps: 

(a)  Preprocessing  by  MDA’s  Pre-Processing  Module 

(b)  Beamforming  using  the  following  CHASP  configuration  options: 


azimuth_focus  No 

azimuth_align  No 

DPCA  Yes 

balance_dpca  Yes 

balance_ratio  0.933 

balance_rad  2.79159 


The  entered  values  arc  examples.  The  value  for  the  phase  (balance_rad)  must  be  k  —  (y 
where  <|)c  defines  the  phase  of  the  steering  vector  G  defined  in  Section  2.2.  Since  the 
DPCA  operation  (channel  subtraction)  is  used,  the  sign  must  be  inverted. 

2.  Creation  of  the  file  with  the  SVs  at  the  rate  of  the  PRF  by: 

(a)  Extracting  a  SV  for  a  known  pulse  of  the  original  data 

(b)  Short  term  orbit  propagation  by  pulse-to-pulse  numerical  integration  using  low-order 
(3)  zonal  harmonics  of  the  Earth  gravitational  potential  (which  could  be  replaced  by  an 
orbital  interpolator) 

3.  Attitude  estimation,  which  includes  the  following  steps: 

(a)  Running  CHASP  in  small  overlapping  tiles  (with  tile  size  of  256  range  bins  and  2048 
azimuth  bins  and  successive  steps  of  128  bins  in  range  and  azimuth)  using  two  iterations 
for: 

i.  Enabling  and  executing  DC  estimation  without  focusing 

ii.  Extracting  preliminary  DC  estimates  from  the  CHASP  header 

iii.  Focusing  with  the  preliminary  DC  estimates  and  performing  DC  estimation  after 
focusing  to  provide  refined  DC  estimates 

(b)  Appending  a  land  mask  flag  to  the  CDC  map  using  georeferencing  information  for  each 
tile  and  using  a  coarse  global  land-water  mask 

(c)  Calculating  Vdown,  vright,  vjor  and  y  defined  in  Section  2.1  for  each  tile 

(d)  Solving  (1)  in  the  least  squares  (LS)  sense 

4.  Information  extraction  from  the  AIS  data,  which  includes: 

(a)  Calculating  image  location  from  the  predicted  ship  geographical  location 
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(b)  Projecting  ship  velocity  given  as  speed  and  heading  onto  the  LOS,  cross-track  direction 
and  along-track  direction 

5.  Extraction  and  analysis  of  each  ship  signature  including: 

(a)  Setting  start  and  stop  range  and  azimuth  sample  for  the  specific  ship  signature 

(b)  Calculating  the  nominal  DC  for  the  range  of  interest  using  the  estimated  attitude  and 
SVs  (applying  (1)) 

(c)  Calculating  the  focusing  parameter  fp/vejj  for  the  range  of  interest  using  the  SVs  (with 
veff  defined  in  [3]) 

(d)  Focusing  the  chip  of  interest  with  the  nominal  parameters  and  saving  the  “before”  image 

(e)  Running  the  tracking  algorithm  iteratively  (four  times)  where  each  iteration  consists  of: 

i.  Enabling  and  executing  the  tracking  algorithm  using  current  processing  values 

ii.  Extracting  the  resulting  DC  and  along-track  speed  correction  from  the  CHASP 
header 

iii.  Setting  the  new  DC  value  and,  optionally,  the  new  along-track  speed  for  the  next 
iteration 

(f)  Saving  the  final  tracking  results  (extracted  Doppler  history  and  power  profiles,  speed 
estimates,  fitting  results) 

(g)  Focusing  with  the  final  DC  value  and,  optionally,  the  along-track  speed  adjustment  and 
saving  the  “after”  image 

(h)  Calculating  target  attributes,  which  includes: 

i.  RCS  and  a°  calculation 

ii.  Region  growing  algorithm 

iii.  Georeferencing 

(i)  Performing  the  map-drift  algorithm  in  steps: 

i.  Zooming  on  the  target  (oversampling  4  times  in  range  and  8  times  in  azimuth)  and 
extracting  two  Doppler  looks 

ii.  Evaluating  look  cross-correlation  and  measuring  azimuth  misregistration 

iii.  Fitting  a  model  to  measured  azimuth  misregistration 

6.  Collection  of  estimation  results  from  all  ships 
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List  of  acronyms 


ASCII 

American  standard  code  for  information  interchange 

AIS 

Automatic  identification  system 

ATI 

Along-track  interferometry 

BW 

Bandwidth 

CHASP 

Chip-based  adaptive  SAR  processor 

CDC 

Clutter  Doppler  centroid 

COASP 

Configurable  airborne  SAR  processor 

DC 

Doppler  centroid 

DOA 

Direction  of  arrival 

DPCA 

Displaced  phase  center  antenna 

DRDC 

Defence  R&D  Canada 

ECEF 

Earth  centred.  Earth  fixed  (system) 

ECI 

Earth  centred  inertial  (system) 

EV 

Eigen-value 

HH 

Horizontal  transmit,  horizontal  receive  (polarization) 

LS 

Least  squares 

LOS 

Line  of  sight 

MDA 

Mac  Donald,  Dettwiler  and  Associates  Ltd.  (Corporation) 

MODEX 

Moving  object  detection  experiment 

PGM 

Portable  gray  map  (format) 

PHP 

Phase  history  product 

PPM 

Portable  pixel  map  (format) 

PRF 

Pulse  repetition  frequency 

R&D 

Research  &  Development 

RCMC 

Range  cell  migration  correction 

RCS 

Radar  cross  section 

RMS 

Root-mean-square 

rmse 

Root-mean-square  error 

SAR 

Synthetic  aperture  radar 

SCR 

Signal-to-clutter  ratio 

SLC 

Single  look  complex 

SNR 

Signal-to-noise  ratio 

vv 

Vertical  transmit,  vertical  receive  (polarization) 

WGS84 

World  geodetic  system  1984 
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